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ABSTRACT 


This  report  records  a  method  of  predicting  the  static  stability 
of  vehicles,  such  as  rough  terrain  forklift  trucks,  on  various  types 
of  slopes  by  computer  analysis. 

Two  basic  methods  are  used  to  obtain  equations  for  determining 
the  critical  slope  for  a  vehicle.  These  are:  (1)  the  action  line 
method,  in  which  the  combined  center  of  gravity  (CCG)  for  the  vehicle 
is  determined,  and  the  critical  slope  obtained  by  finding  the  sides lope 
upon  which  the  vehicle  must  be  resting  so  that  the  CCG  is  directly 
over  the  action  line  formed  by  the  two  downhill  points  of  support  of 
the  vehicle,  and  (2)  the  wheel  load  method,  in  which  the  loads  on  the 
four  tires  are  examined  under  all  possible  sideslope  conditions  to 
determine  the  minimum  slops  for  which  the  vehicle  will  be  in  an 
unstable  condition. 

The  report  includes  *  computer  program  using  the  equations  derived 
from  the  two  methods  for  determining  critical  slopes.  This  program 
allows  the  vehicle  parameters  such  as  type  of  steer,  suspension,  frame, 
weights,  and  dimensions,  to  be  varied,  and  for  each  set  of  parameters 
provides  the  maximum  slope  on  which  the  vehicle  can  rest  in  a  stable 
condition.  The  program  also  shows  the  orientation  of  the  vehicle 
corresponding  to  this  critical  slope. 

The  computer  program  follows  the  wheel  load  method  for  vehicles 
with  midrange  oscillation;  i.e. ,  vehicles  in  which  the  front  part  can 
rotate  relative  to  the  rear  part  about  a  longitudinal  axis,  with  the 
oscillation  joint  located  somewhere  between  the  front  and  rear  axles. 
For  all  other  types  of  vehicles,  it  uses  the  action  line  method. 


cOmfuiisk  ANALYSIS  OF  FORKLIFT  TRUCK  STABILITY  WHEN 
OPERATING  ON  SIDE  SLOPES  UNDER  NEAR  STATIC  CONDITIONS 


1 .  Introduction 


The  objective  of  this  program  is  to  determine  the  maximum  side 
slope  and  orientation  on  which  a  theoretical  or  actual  forklift  truck 
can  safely  operate.  A  commercially  available  articulated  forklift 
truck  was  evaluated  to  determine  its  potential  for  army  materials 
handling  use  in  the  field.  While  investigating  the  stability  of  this 
vehicle,  the  theoretical  maximum  sideslope  was  desired.  This  presented 
no  problem  when  the  vehicle  was  operating  straight  across  the  sideslope 
parallel  to  the  bottom  of  the  slope.  The  CCG  could  be  calculated  from 
the  CCG  of  the  vehicle  and  load.  In  addition,  one  could  readily 
determine  the  angle  or  percent  of  slope  required  to  move  the  CCG  to 
the  action  line  formed  by  the  contact  points  of  the  two  downhill  tires. 
The  problem  occurred  when  the  vehicle  turned  up  or  down  slope  due  to 
the  geometry  of  it’s  articulated  frame  and  articulated  steering.  A 
small  degree  of  steering  would  shift  both  the  center  of  gravity  (CG) 
of  the  front  mass  of  the  vehicle  and  the  load  considerably,  in  an 
x,  y  and  z  direction  as  compared  to  a  vehicle  with  Ackermann  steering 
(common  automobile-type  steering)  and  a  fixed  frame. 

The  individual  CG's  of  the  separate  masses  were  difficult  to 
determine  as  the  vehicle  moved  in  a  circle  on  the  sideslope  increasing 
its  percentage. 

It  was  iecided  that  due  to  the  number  of  variables  involved,  a 
computer  program  could  be  written  to  include  not  only  most  of  the 
constants  and  variables  pertaining  to  the  vehicle,  but  would  allow 
these  constants  and  variables  to  be  changed  as  the  conditions  required, 
and  would  locate  the  CCG  relative  to  the  action  line  of  tne  vehicle. 

Further  investigation  indicated  that  this  program  should  include 
vehicles  with  rigid  suspension  and  different  types  of  frames  and  steer¬ 
ing.  There  are  two  analytical  methods  of  determining  vehicle  stability 
on  side  slopes.  One  method,  which  we  shall  call  Case  I,  is  to  place 
the  vehicle  theoretically  on  a  side  slope  and  rotate  the  vehicle  from 
0°  to  360°  and  increase  the  slope  until  the  CCG  falls  on  the  action 
line  of  the  downhill  tires. 
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As  the  CCG  approaches  the  action  line,  the  vehicle  becomes  more  unstable 
and  finally  turns  over  laterally  when  the  CCG  passes  beyond  its  limits. 
This  determines  both  the  maximum  side  slope  on  which  vehicle  can  operate 
and  its  orientation. 

The  second  method,  Case  II,  is  to  place  the  vehicle  on  a 
horizontal  plane  and  tilt  tne  plane  in  a  given  directio.-  until  one  or 
more  of  the  two  uphill  tires  has  a  zero  wheel  loading.  The  direction 
of  down  slope  is  changed  1°  at  a  time  and  the  process  repeated  until 
critical  slopes  are  obtained  for  downhill  directions  of  0°  through  360°. 
The  critical  slope  and  orientation  of  the  vehicle  will  be  obtained  by 
selecting  the  minimum  of  the  critical  slopes  with  its  corresponding 
down  slope  direction.  This  will  be  the  maximum  side  slope  and  orientation 
on  which  the  vehicle  can  safely  operate. 

These  two  cases  or  methods  are  followed  in  this  report.  Case  I  is 
used  for  vehicles  with  articulated  steer,  articulated  frame,  and  rear 
axle  oscillation;  Ackermann  steer,  straight  frame  and  rear  axle  oscil¬ 
lation;  and  any  vehicle  in  which  no  oscillation  is  possible.  Case  II 
is  used  for  vehicles  with  articulated  steer,  articulated  frame,  and 
midrange  oscillation;  and  rear  axle  steer  (wagon  or  Ackermann),  straight 
frame  amd  mid-range  oscillation.  This  report  (Phase  T)  will  handle  only 
the  static  conditions  where  the  vehicle  is  stationary  or  barely  moving 
on  the  side  slope.  A  second  report  (Phase  II),  planned  for  a  later 
date,  will  cover  the  dynamic  conditions  where  vehicle  speed  (centrifugal 
force)  wheel  torque,  and  soil  conditions  will  be  considered. 

2 .  Case  I .  The  Action  Line  Method  of  Determining  Critical  Slope  and 
Orientation 


a.  The  Equations  for  Critical  Slope  and  Orientation 

In  determining  the  minimum  degree  of  slope  on  which  a  vehicle  is 
unstable,  it  is  assumed  that  the  only  forces  acting  on  the  vehicle  are 
those  caused  by  the  weight  of  the  vehicle  and  the  load. 

In  using  the  action  line  method  of  accomplishing  the  above,  we 
determine  the  minimum  slope  (and  corresponding  direction  of  downslope) 
for  which  a  given  CCG  is  directly  over  a  given  action  line,  the  action 
line  being  a  line  through  two  points  of  support  of  the  mass. 
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Suppose  an  object  with  a  center  of  gravity  (CG)  at  the  point 
C  =  <xc»yc»2c>  is  supporter.’  -at  three  or  more  points  P^  =  (  xi ,  Yi ,  zi ) , 
**2  ~  ^x2*y2»z2^»  •  *  *  *as  sl*own  in  Figure  1. 


Z 


Figure  1.  The  Suspension  Triangle 

Assume  that  the  x,y  plane  is  horizon ted.  and  that  the  points  C,  Plf  P2,.., 
are  fixed  to  the  given  axis  system,  but  none  of  the  points  need  lie  in  the 
horizontal  plane.  (A  top  view  is  shown  in  Figure  2. )  If  the  whole  system 
is  tilted,  so  that  the  x,y  plane  is  no  longer  horizontal  but  has  a  downslope 
direction  as  indicated  in  Figure  2,  the  point  C  will  move  in  a  vertical 
plane  p  parallel  to  the  downslope  direction.  If  the  system  is  tilted  far 
enough  in  the  given  direction,  the  CG  will  eventually  be  directly  over 
the  point,  P,  which  is  the  intersection  of  the  vertical  plane  p  with  the 
line  P^. 
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Figure  2.  Top  View  of  Suspension  Triangle 

It  is  evident  that  the  angle  through  which  the  system  must  be 
tilted  in  order  that  the  point  C  be  directly  over  the  point  P,  is  the 
angle  0  formed  by  the  segment  PC  and  a  vertical  ray  from  P,  as  shown 
in  Figure  3.  Hence,  if  the  system  is  tilted  through  an  angle  0  in  the 
given  direction,  an  unstable  condition  will  exist,  since  the  horizontal 
projection  of  the  CG  will  no  longer  lie  in  the  interior  of  the  horizontal 
projection  of  the  support  triangle.  (In  the  figures  shown,  it  is  a 
triangle,  but  it  may  be  a  quadrilateral  or  some  other  configuration, 
depending  on  the  type  of  suspension  of  the  vehicle. ) 

In  the  following  discussion,  we  shall  use  the  notation: 


pl  = 

yx,  z !>» 

P2  = 

(X2’  y2’  V* 

C  = 

(xc»  yc> 

Ax  = 

X2  "  Xl* 

il 

>> 

<1 

*2  -  *1. 

Az  = 

Z2"  V 

4 


C  (c.c.q.) 


VTP» 


Vertical  ray 


•Vertical  plane  p 


Figure  3.  The  Critical  Slope  Angle,  0 

We  will  use  brackets  ([  J)  to  denote  vectors;  i.e.,  A  = 

£alt  a2,  a^Q  denotes  tne  vector  with  x,  y,  and  z  components 

ax,  a2,  and  ag,  respectively. 

In  Figure  3,  the  point  P  =  (x,  y,  z)  is  on  the  line  containing  and  P, 
Any  such  point  P  can  be  written  in  the  form: 

rx  =  xx  +  t(Ax), 

(1)  '  y  =  yx  +  t(Ay), 

z  =  z  +  t(Az). 


The  value  of  t  determines  the  position  of  P  -  (x,  y,  z)  relative  to 

P  and  P  ;  i.e.,  if  t<0,  P  is  on  the  "dotted”  portion  of  the  line 
1  2 

(Fig.  4);  if  t  =  0,  P  =  Pxj 


Figure  4.  The  Action  Line 

if  0<t<l,  P  is  between  Pj  and  P2;  if  t  =  1,  P  =  P2;  and  finally  if 
t»l,  P  is  in  the  "dashed"  portion  of  the  line. 


We  wish  to  find  an  analytical  expression  for  the  angle  0  in  Figure  3. 
We  recall  from  vector  analysis  that  the  angle  0  formed  by  two  vectors 


A  and  B  is  given  by: 

(2)  cos  0 


— »  — > 


IA|  IB  I 


where  A*B  is  the  dot  product,  and  the  absolute  value  symbols  denote 
vector  magnitudes.  Hence,  the  angle  0  in  Figure  3  is  given  by; 


(3)  cos  0 


(PC)  *  (Z) 

—  -  '  ■> 

|  PC|  \z\ 


where  PC  =  |x  ~  x>  yc“  y»  zc  -  zjis  the  vector  from  P  to  C,  and 
Z  =  [o,  0,  ljis  the  unit  vertical  vector.  Using  the  fact  thatIZ  1=  1, 

equation  (3)  becomes: 
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u? 


(4) 


cos  0 


<*c  -  z> 


/ 


(xc  -  x)2  I  (yc  -  y)2  i-  (zc  -  z)‘ 


Substituting  equations  (1)  into  (4),  we  have: 
(5)  cos  0  =  _ 


z  -  z.  -  tAZ 
£  *  — 


^  <V  X1- 


2  2  2 
tAx)  +  (yc-y1~tAy)  +  ( z^-z-,  -tAz) 


~‘c  1 


Equation  (5)  gives  the  angle  0  for  a  given  value  of  t;  i.e.,  for  a 
given  point  P  on  the  line  containing  Pj  and  p£.  We  wish,  however,  to 
find  the  minimum  angle  0  as  P  ranges  over  the  given  line.  We  can 
accomplish  this  by  maximizing  the  function  of  t: 


f(t)  -  cos  0  - 


Zc  -  Zj  -  tAZ 


y/(xc-  x  tAx)2  7 


2  2 
*•  (yc-yi'tAy)  +  (Zp-zj-tAz) 


The  function  f(t)  is  of  the  form: 


(6)  f(t) 


a  +  bt 


y/ 


c  +  dt  -  et* 


where  a  =  z  -  z, 

c  1 

b  =  -Az 

c  =  (xc  -  xx)2  +  (yc  -  yi)2  +  (zc  -  zx)2 
d  -  -2[axCxc  -  x1>  +  Ay(yc  -  +  Az(zc  -  zx)J 

e  =  (ax)2  +  (Ay)2  *■  (Az)2. 

Using  the  usual  calculus  methods  to  maximize  f(t),  we  set  f’(t)  = 
0  and  solve  for  t: 


(7)  f’(t)  =  _ b _  -  (a-fot)  .  (d+2et) 


V 


c+dt+et2 


(c'dti ct2) 


0, 


7 


(8) 


i 


t 


Hence,  if  we  let  tQ 


ad  -  2b c 
bd  -  2ac 

ad  _ 2  be 

bd  -  2ac 


then  the  relation 


(9) 


cos  0  =  a  +  bt, 


o_ 


y/ c  +  dtQ  +  etQ2 


gives  the  minimum  angle  0  as  P  ranges  over  the  whole  line  containing 

and  P2#  Note  that  we  are  not  interested  in  the  whole  line,  but  only  in 
the  segment  PjP2-  That  is,  at  this  time  we  are  interested  only  in  downslope 
directions  for  which  the  corresponding  vertical  plane  p,  shown  in  Figure  3, 
intersects  the  line  segment  PjP2-  If  t(j>l,  which  corresponds  to  a  point 
outside  the  segment  P^^.the  minimum  0  for  the  range  in  which  we  are 
interested  will  occur  at  t=l.  Similarly,  if  tQ^b  the  minimum  0  for  the 
range  in  which  we  are  interested  will  occur  at  t  =  0.  The  previous  state¬ 
ments  follow  from  the  fact  that  the  function  0  =  0  (t)  is  monotonically 
increasing  as  t  proceeds  in  either  direction  from  tQ. 


Ckice  we  have  determined  the  value  t  =  tg  corresponding  to  the  minimum 
slope  angle,  we  can  easily  find  the  corresponding  downslope  direction  by 
finding  the  horizontal  projection  of  the  vector  C&,  which  lies  in  the 
vertical  plane  p  as  can  be  seen  from  Figures  2  and  3.  For  t  =  tQ  this 
horizontal  projection  is  the  vector: 

v=  fcc-Xc,  y-yc,  0j=  [xi+t0Ax-xc,  yi+tt/y-y,  o) . 

If  0  is  the  angle  formed  by  the  X-axis  and  the  horizontal  vector  V,  then 
0  can  be  determined  from  the  relation: 


(10)  tan  0  =  yi  +  t0  liy  -  yr 

Xj  +  t0Ax  -  xc 


By  using  the  equations  obtained  in  the  preceding  discussion,  one  can 
easily  find  the  critical  slope  and  downslope  direction  relative  to  a  given 
OG  and  action  line.  Usually  one  must  consider  several  action  lines  when 
studying  the  static  stability  properties  of  a  given  object,  the  action 
lines  being  determined  by  the  points  of  support  for  the  object  in  question. 
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The  overall  critical  slope  would  be  the  minimum  of  the  critical  slopes 
relative  to  the  various  action  lines,  and  the  overall  critical  down 
slope  direction  would  be  the  down  slope  direction  corresponding  to  the 
overall  critical  slope. 

b.  The  Action  Line  for  Rubber-Tired  Vehicles 


A  vehicle  with  low  pressure  rubber  tires  is  not  supported 
by  contact  "points'*  but  rather  by  contact  "areas,"  namely,  the  tire 
"footprints."  If  we  are  to  use  the  action  line  method  of  determining 
critical  slope  for  such  a  vehicle,  we  must  make  a  reasonable  assumption 
of  what  point  in  the  footprint  of  a  tire  to  take  for  the  end  point  of  an 
action  line  (Fig.  5). 


Pi 


Tire  Footprints 


Figure  5.  Tire  "Footprints" 

One  method  might  be  to  take  the  center  of  the  footprint  for  an  end  point 
of  the  action  line.  However,  we  must  remember  that  the  true  action  line 
is  the  line  about  which  the  vehicle  will  rotate  once  the  COG  passes  over 
this  line.  Consider  a  vehicle  on  a  sideslope  and  suppose  that  the  vehicle 
is  in  an  unstable  condition;  that  is,  at  this  moment  the  vehicle  has  not 
tipped  over,  but  if  the  slope  were  increased  by  a  small  amount  (in  the 
existing  downslope  direction),  turnover  would  occur.  In  this  case,  the 
vehicle  would  be  balanced  on  the  downhill  tires  with  zero  load  on  the  uphill 
tires.  The  situation  would  be  similar  to  that  depicted  in  Figure  6. 

If  we  take  moments  about  a  line  L  directly  under  the  COG  and  passing 
through  the  tire  footprints  (Fig.  7),  the  resultant  moment  about  this  line 
must  be  zero  since  the  vehicle  is  assumed  to  be  in  equilibrium. 
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L,et 


A  and  is  denote  "the  respective  project! who  of  ths  «u«sl  centers 


onto  the  plane  of  the  slope.  Then,  for  our  assumed  equilibrium 
condition  the  line  segment  AB  cannot  lie  directly  under  the 


4*  CG 


Figure  6.  Force  Distribution  on  Tire  Profile 

COG,  for  if  AB  did  lie  directly  under  the  COG,  the  moments  tending  to 
cause  turnover  about  the  line  AB  would  be  greater  than  the  moments 
acting  in  the  opposite  direction.  This  would  happen  since  the 
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compression  on  the  uphill  side  of  the  tire  (Fig.  6)  causes  the  pressures 
uphill  from  AB  to  be  much  greater  than  the  pressures  downhill  from  AB. 

The  above  discussion  indicates  that  if  the  sideslope  were  increased, 
turnover  would  occur  about  a  line  uphill  from  the  segment  AB.  Exactly 
how  far  uphill  this  action  line  would  be  located  depends  on  the  elastic 
properties  of  the  tire  as  well  as  the  weight  of  the  vehicle.  Lack  of 
complete  information  on  the  elastic  properties  of  tires  makes  it  difficult 
to  predict  the  true  position  of  the  action  line,  but  it  seems  reasonable 
that  the  end  points  of  the  action  line  should  be  taken  somewhere  between 
the  center  and  the  inside  edge  of  the  tread. 


4 


Downs lope 


Figure  7.  Actual  Action  Line 

3.  Case  II.  Wheel  Load  Method  of  Determining  Critical  Slope  and 
Orientation 


a.  General  Discussion 


In  using  the  wheel  load  method  for  studying  the  stability  of  a 
vehicle,  we  determine  the  loads  on  each  of  the  four  tires  of  a  vehicle 
with  a  given  orientation  on  a  given  side  slope,  and  decided  whether  or 
not  the  vehicle  is  in  equilibrium.  If  we  test  the  vehicle  for  equilibrium 
at  all  possible  side  slopes  and  orientations,  we  can  then  choose  the 
minimum  side  slope  and  corresponding  orientation  for  which  the 
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vehicle  is  not  in  equilibrium.  This  will  give  the  critical  slope  and 
orientation  for  the  vehicle  under  consideration. 

Usually  the  wheel  loads  (in  a  condition  of  equilibrium)  are 
determined  by  solving  the  equilibrium  equations  for  the  unknown  forces 
on  the  tires.  For  a  given  side  slope  and  orientation,  we  will  assume 
that  the  vehicle  is  in  equilibrium  and  compute  the  wheel  loads  from 
the  equilibrium  equations.  The  forces  on  the  tires  obtained  will  be 
the  forces  necessary  for  equilibrium  to  exist.  If,  in  reality,  the 
vehicle  is  not  in  equilibrium,  one  of  the  computed  forces  will  be  acting 
in  a  direction  physically  impossible  for  our  system;  i.e. ,  our  results 
will  tell  us  that  the  system  will  be  in  equilibrium  provided  we  apply 
am  additional  downward  force  to  one  of  the  uphill  tires,  which  is 
equivalent  to  saying  the  system  is  not  in  equilibrium  as  it  stands. 

In  using  the  wheel  load  method  we  shall  consider  vehicles  in  which 
midrange  oscillation  is  possible;  i.e.,  at  some  point  between  the  front 
and  rear  axles  the  front  part  of  the  vehicle  cam  rotate,  relative  to  the 
rear  part  of  the  vehicle,  about  some  fixed  axis .  Assume  the  vehicle  is 
composed  of  two  rigid  bodies  joined  at  one  point  J  so  that  the  two  bodies 
cam  rotate  relative  to  one  amother  about  some  fixed  axis  L  passing 
through  the  common  point  (Fig.  8); 


Front 


oxit 


Figure  8.  Suspension  Triangles 
for  Vehicle  with  Midrange  Oscillation 


12 


Since  we  are  concerned  here  with  only  static  stability  where  the 
only  external  forces  are  those  caused  by  the  weight  of  the  vehicle,  we  will 
assume  that  the  resultant  reaction  forces  acting  on  each  of  the  four  tires 
are  vertical.  In  the  following  analysis,  we  will  assume  that  the  coordinate 
axes  are  chosen  in  some  convenient  manner  so  that  the  x,y  plane  is  horizontal. 
Assume  that  the  vehicle  is  resting  on  some  plane,  with  the  tire  loads 
centered  at  the  points  R.  (x.^,  yi ,  z..),  i  3,  2,  3,  4,  where  Rj  and  R0 
refer  to  the  front  tix*es  and  R3  and  H4  refer  to  the  rear  tires.  Our 
notation  will  be  as  follows: 

w^  “  weight  of  the  front  part  of  the  vehicle,  including 

load, 


RCl  =  (xCi>  yc  ,  zCl>  =  center  of  gravity  of  the  front  part  of 
the  vehicle, 

W2  =  weight  of  the  rear  part  of  the  vehicle, 

RC2  =  (xc2>  yC2«  zc2^  =  center  of  gravity  of  the  rear  part 
of  the  vehicle, 


-1-  w. 


total  weight  of  the  vehicle. 


Rc  =  Cxc,  yc,  zc)  =  combined  center  of  gravity  of  the  vehicle 
where  xc  Xci  (  wi/w)  +  xc^  (w2/w),  yc  yCi  (wx/w)  :  yc^  (w2/w),  and 

zc  ~  zCi  <wi/w)  +  zC2  (W2/W), 

Pi  =  [0,  0,  fJ]  =  the  veexor  representing  the  vertical  reaction 
force  acting  on  the  tire  at  the  point  Rj;  i  =  1,  2,  3,  4, 


Wi  =  £0,  0,  -wj|  =  the  vector  representing  the  weights  of  front 

and  rear  parts  of  the  vehicle,  acting  at  the  points  R_  respectively; 

ci 

i  =  1,  2, 


W  =  |_0,  0,  -wj  =  the  vector  representing  the  total  weight  of 
the  vehicle,  acting  at  the  point  Rc< 

— b 

If  A  denotes  the  point  <a^,  a2,  a2),  then  A  will  denote  the 
vector  j  a  ,  a  a  “j  representing  the  directed  segment  from  the  origin 
to  the  point  A.  We  i*eca3 1  from  the  Principles  of  Mechanicsv  - 1  that  if 


F  is  a  force  acting  at  a  point  B,  then  the  vector  moment  of  F  about  a 
point  A  is  given  by: 

M  =  ABXF 

where  AB  is  the  vector  from  A  to  B  and  the  symbol  :,X’'  denotes  the  vector 
cross  product  (Fig.  9).  We  recall  also  that  if  L  is  a  line  containing 
A  and  U  is  a  unit  vector  parallel  to  L,  then  the  scalar  moment  of  "p* 
about  the  line  L  (with  the 


Figure  9.  Moment  about  a  Point 
positive  direction  determined  by  U)  is  given  by  (Fig.  10). 


M  =  U  •  (AB  X  F) 


Figure  10.  Moment  about  a  Line 
b.  Wheel-Load  Equations 

We  shall  now  obtain  the  equations  for  the  wheel  loads  , 
i  =  1,  2,  3,  4.  Assume  that  our  vehicle  is  resting  on  a  slightly 
tilted  plane,  in  a  condition  o'  equilibrium. 
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/  o  \ 

A  necessary  condition'- “  'or  the  system  to  be  in  equilibrium  is  that: 

( 1 )  "f  =  ~0  and 

(2)  = 


where  F  is  the  sum  of  the  applied  forces  acting  on  the  system 

and  G  is  the  sum  of  the  applied  moments  about  some  fixed  point,  say  the 

origin.  For  our  system,  equations  (1)  and  (2)  are  respectively: 

o)  ~f  =  r  'Fi  +  "W  =0, 

i  1 
4 

(4)  G  ~  Z  \l±  X  Fi  -r  Rc  x  W  -  0. 

i  1 


Upon  examining  the  individual  components  of  the  V'ectors  in 
equations  (3)  and  (4),  we  find  that  (3)  reduces  to  the  scalar  equation 

(5)  22  Fi  =  w, 

i=l 


and  (4)  reduces  to  the  two  scalar  equations 


(6)  < 


£  xi  Fj  =  w  xc, 
i -1 


4 

s  n  Fi 

i  1 


w  y 


c* 


Our  system  of  equations  for  the  determination  of  the  F  is: 

i 


(7)  -A 


F  F  +  F  +  F  -"w 

1  2  3  4 


x-  F,  +  x„  Fn  +  x0  F„  +  x.  F.  =  w  x 


J  1  1  2  2 


y  F  +  y  F 
*11  *2  2 


3  3  4  *4 


+  y„  F  +  y  F  =  w  y  < 

J  3  4  4  C 
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We  can  assume  that  the  three  equations  in  (7)  are  independent  since  a 
3X3  determinant  from  the  coefficient  matrix 

"l  1  1  1 

X1  x2  X3  X4 

yi  y2  V3  V4 

is  zero  only  if  three  of  the  points  R^,  R2,  R3,  R^  are  collinear 
which  certainly  will  not  be  the  case  for  any  vehicle  we  wish  to  consider. 
Any  three  of  the  can  be  obtained  in  terms  of  the  fourth  F^  from  (7); 
however,  we  cannot  get  a  unique  solution  from  the  system  of  three  equations 
in  four  unknowns.  Therefore,  we  have  not  supplied  enough  information  to 
uniquely  determine  the  F^,  and  we  must  look  for  a  fourth  equation  in  the 

F^.  The  fourth  equation  comes  from  the  fact  tnat  the  vehicle  is  composed 
of  two  rigid  bodies  which  can  rotate  relative  to  one  anotier  about  a 
fixed  axis  L.  If  the  vehicle  is  to  be  in  equilibrium,  then  the  resultant 
scalar  moment  of  the  front  part  of  the  vehicle  about  the  axis  L  must  be 
zero.  If  J  =  (xj,  yj,  zj)  is  the  common  point  in  Figure  8,  and 
U  =Julf  u2,  u^jjis  a  unit  vector  parallel  to  L  and  directed  toward  the 
front  of  the  vehicle,  then  the  resultant  scalar  moment  01  the  front  part 
of  the  vehicle  about  L  (with  positive  direction  determined  by  U)  is  given 
by 

(8)  m  =  U» ( JRj  X  F^ )  +  x  +  "U'C^c  X~W^) 


«1  u2  u3 

U1  u2  u3 

U1  u2  u3 

xi~xj  yi-yj  zi~zj 

+ 

x2-xj  y2-yj  z2_zj 

+ 

xCi~xj  yCjL-yj  zCi"zj 

0  0  Fx 

0  0  f2 

0  0  -Wj 

U1  u2 

U1  u2 

U1  u2 

=  F1 

xi~xj  yi-yj 

+  F2 

x2~xj  y2~yj 

—  w- 
1 

xC;l-xj  yCl~yj 

=  Fxdx  +  F2d2  -wxd 


,  where 


u 

u 

1 

2 

dl = 

I  xrxj  yryj 

f 

d2 

U1  u2 


u 

*> 


y2~-vj 


i 

! ,  and 
I 


Setting  the  scalar  moment  m  in  (8)  equal  to  zero,  we  obtain  the  equation 


(9) 


dlFl  + 


a2F2 


w. 


Hence,  (7)  and  (9)  give  the  following  system  of  four  equations 
in  four  unknowns : 


dlFl  - 

d2F2 

w,  d 

1  C1 

F1  ‘ 

r2  h 

F3 

+  F4  " 

w 

(10)  S 

X1F1 

X2F2  h 

X3F3 

•'  x  .F„  “ 
4  4 

W  X 

c 

yiFl  + 

y2F2  + 

y3F3 

+  y  F  = 
y4  4 

w  y 

c. 

Note  that  equations  (10)  force  the  resultant  scalar  moment  of 
the  rear  pare  of  the  vehicle  about  L  to  be  zero  also,  since  the  forward 
moments  about  L  are  zero  and  since  equations  (1)  and  (2)  require  that 
the  cum  of  all  the  moments  about  L  be  zero. 


Equations  (10)  can  easily  be  solved  for  the  Fj,  by  means  of 
determinants.  The  solutions  to  (10)  are  as  follows: 


(11)  Fj  =  D^/D  for  i=l,2,3,4,  where 

D  --  determinant  of  coefficients 

=  dl  l/W  ~  (x4'x2}  (y3"y2> 

-  d0  f(x0-x,)  (y4-yx)  -  (x„-xn)  Cy^-y,)' 


3  1 


4  1 


3 


D,  =  w,  d 
1  1  c. 


-  w  d. 


[<x3"x2> 

[<W  ‘W 


(xrx2)  (y3-y2)] 
‘V'c1  (y3‘yc>], 


1? 


d2  =  *  di  DvV  ‘W  -  ‘W <  vv] 

'  V^C'VV  <y4'yl>  ■  (W  <y3_yl)]' 

°3  =  W  dl  [(xc"x2>  (y4~y2}  "  (x4~x2}  (yc~y2) 

~  W  d2  QvV  (y4"yl>  "  (X4"X1)  (Vyl>; 

+  wldc1[<X2~Xl)  (y4~yl)  "  (X4"X1)  (y2-yl}]- 

D4  =  w  ^  [(x3-x2)  (yc-  y2)  -  (xc-x2)  (y3-y2)] 

-  w  d2  [(x3-xx)  (yc-  yx)  -  (xc-Xl)  (y^)] 

-  w1dcJ(x2-x1)  (y3-  yx)  -  (x3-xx)  (y2-yi)]. 

The  unit  vector'd  =  ,  u2,  uj,  in  (8),  which  is  parallel 

to  the  oscillation  axis  L,  can  be  determined  if  we  know  two  points 
P  and  Q  on  the^line  L.  If  P  =  <xp,  yp,  Zp)  and  Q  =  (x^,  y^,  zq) 
are  on  L  and  PQ  is  directed  toward  the  front  of  the  vehicle,  tnen 
we  can  let  U  be  the  unit  vector: 

~U  =  (1/jPQj)  PQ 

|^Q_xp*  yQ-yP '  ZQ_ZpJ* 


After  examining  equations  (8),  (9),  and  (11)  we  see  that  the  factor 

(1/|PQ|)  in  U  will  appear  as  a  factor  in  the  numerator  and  denominator 

in  (11)  and  hence  will  cancel  out.  Therefore,  to  simplify  calculations, 

we  can  take  U  =  PQ  to  determine  d,  ,  d_,  and  d 

12  c. 
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4.  Description  of  Computer  Program 

A  computer  program  which  calculates  critical  slopes  and 
orientations  for  various  types  of  rigid-suspension  vehicles,  uses  the 
action  line  method  (Case  I)  for  the  following  types  of  vehicles: 

(a)  articulated  steer,  articulated  frame,  and  rear  axle  oscillation, 

(b)  Ackermann  steer,  straight  frame,  and  rear  axle  oscillation,  and 

(c)  any  rigid-suspension  vehicle  on  which  no  oscillation  is  possible. 

The  program  uses  the  wheel  load  method  (Case  II)  for  any  rigid- 
suspension  vehicle  having  midrange  oscillation  (as  described  in  Section  3). 

In  case  oscillation  can  occur,  the  program  computes  critical  slopes 
for  oscillation  and  turnover.  A  general  description  of  the  computer 
program  is  given  in  the  following  paragraphs.  If  more  details  are 
desired,  a  listing  of  the  Fortran  program  is  given  in  Appendix  IV. 

A  general  flow  diagram  is  shown  on  page  30  in  this  Section  (Fig.  13). 

TI»e  basic  notation  used  in  the  computer  program  is  shown  in 
Figure 


Figure  11.  Standard  Position  for  Input  Data 

The  axic  system  is  assumed  to  be  oriented  so  that  the  vehicle  is 
resting  on  the  x,y  plane  with  tne  x-axis  parallel  to  the  longit  final 
axis  of  the  vehicle,  and  the  origin  is  at  the  projection  of  the  center 
of  the  rear  axle  onto  the  x,y  plane. 
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The  variable  names  as  used  in  the  program,  and  their  meanings  are 
as  follows: 

FL,  FR,  RR,  RL:  Centers  of  the  bearing  surfaces  of  the  four  tires: 
front  left,  front  right,  rear  right,  and  rear  left,  respectively. 

AJT;  Center  of  the  articulation  joint. 

OJT:  Center  of  the  oscillation  joint. 

FAX:  Center  of  the  front  axle. 

CF:  CG  of  the  front  mass  of  the  vehicle. 

WF:  Weight  of  the  front  part  of  the  vehicle. 

CLD:  CG  of  the  load. 

WLD:  Weight  of  the  load. 

CXF:  OG  of  any  extra  mass  added  to  the  front  part  of  the 
vehicle. 

WXF:  Weight  of  any  extra  mass  added  to  the  front  part  of  the 
vehicle. 

CR:  OG  of  the  rear  mass  of  the  vehicle  (minus  the  rear  axle). 

WR:  Weight  of  the  rear  part  of  the  vehicle  (minus  the  rear 

axle). 

CRAX;  Center  of  the  rear  axle  (also  the  CG  of  the  rear  axle). 

WRAX:  Weight  of  the  rear  axle. 

CXR:  CG  of  any  extra  mass  added  to  the  rear  part  of  the  vehicle. 

WXR:  Weight  of  any  extra  mass  added  to  the  rear  part  of  the 
vehicle. 

ALPHA:  Steer  angle  (always  assumed  to  be  such  so  that  the 
vehicle  will  turn  left  if  ALPHA  is  positive). 
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OSC:  The  oscillation  angle;  i.e.,  the  angle  through  which  the 
front  (or  rear)  part  of  the  vehicle  can  rotate  from  the  neutral  position 
before  hitting  the  mechanical  stop. 

SHIFT:  Distance  from  the  center  of  the  tire  tread  (toward  the 
inside  edge  of  the  tread)  so  that  the  endpoint  of  the  action  line  will 
be  located. 

NSTEER:  Equals  1  if  the  vehicle  has  articulated  steer;  2  if 
rear  wagon  wheel;  and  3  if  Ackermann. 

NOSC:  Equals  1  if  no  oscillation  is  possible;  2  if  rear  axle 
oscillation  is  possible;  and  3  if  midrange  oscillation  is  possible. 

Weights  and  distances  can  be  expressed  in  any  convenient  units 
as  long  as  the  same  units  are  used  throughout.  However,  it  is  suggested 
that  weights  be  expressed  in  lbs  since  the  output  column  headings 
are  labeled  "lbs*’  for  Case  II  vehicles. 

The  symbols  defined  above  which  represent  locations  (in  space) 
are  actually  arrays .each  consisting  of  the  three  coordinates  of  the 
location.  Details  on  the  formats  and  arrangement  of  the  input  data 
cards  are  given  jn  Section  4.d. 

During  execution  of  the  computer  program,  the  computer  first 
reads  in  the  data  for  the  vehicle  in  standard  position  as  shown  in 
Figure  11.  The  coordinates  are  then  recomputed  as  necessary  depending 
on  the  type  of  steer  and  angle  of  steer.  The  computer  then  follows 
one  of  two  possible  sequences  of  events  depending  on  whether  or  rot 
midrange  oscillation  is  possible. 

a.  Case  I  Vehicles 


This  part  of  the  program  pertains  to  vehicles  which  do 
not  have  midrange  oscillation. 

(1)  For  a  vehicle  in  which  no  oscillation  is  possible, 
the  computer  first  calculates  the  combined  CG  of  the  whole  vehicle, 
using  the  usual  moment  equations  which  can  be  obtained  from  any  mechanics 
testbook.  The  computer  then  calculates  the  critical  slopes  relative  to 
each  of  three  action  lines;  that  formed  by  the  left  two  tires;  by  the 
right  two  tires,  and  by  the  front  two  tires.  It  is  left  to  the  person 
using  the  program  to  decide  which  of  the  three  slopes  is  the  minimum. 
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(2)  For  a  vehicle  with  rear  axle  oscillation,  the 

1  9 

configuration  will  be  as  shown  in  Figure  12  (shown  for  a  vehicle  with 
articulated  steering).  In  this  case,  the  combined  CG  of  the  whole 


Figure  12.  Suspension  Triangle  for  Rear  Axle  Oscillation 

vehicle  is  computed,  and  critical  slope  and  orientation  for  turnover 
about  the  left  tires  is  obtained  as  though  the  vehicle  were  of  the  non¬ 
oscillating  type  (mainly  for  comparison  purposes  since  oscillation  will 
probably  occur  before  turnover). 

Next  the  combined  CG  of  the  vehicle  minus  the  rear  axle  is 
computed,  and  critical  slope  and  orientation  are  obtained  using  this  CG 
and  the  action  line  formed  by  the  center  of  the  rear  axle  and  the  left 
front  tire.  This  gives  the  critical  slope  S1  for  osrci'  lation  to  the  left; 
i.e.,  the  minimum  slope  and  corresponding  downslope  d.rection  for  which 
the  CG  for  the  vehicle  minus  the  rear  axle  will  be  over  the  action  line 
RB  in  Figure  12. 

It  is  then  assumed  that  complete  oscillation  to  the  left  occurs; 
that  is,  the  vehicle  oscillates  to  the  left  until  it  hits  the  mechanical 
stop.  The  question  now  is:  In  oscillating  to  the  left,  has  the  combined 
CG  for  the  whole  vehicle  been  thrown  over  the  action  line  AB?  This 
question  is  answered  by  recomputing  the  combined  CG  for  the  whole  vehicle, 
assuming  the  vehicle  is  completely  oscillated  to  the  left,  and  using  this 
CG  to  obtain  the  critical  slope  S2  relative  to  the  action  line  AB. 
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Then ,  if  S^j  is  greater  than  ,  turnover  could  not  have  occurred  when 
the  vehicle  oscillated,  for  the  critical  slope  for  an  oscillated  vehicle 
is,  in  this  case,  greater  than  the  slope  at  which  oscillation  occurs.  On 
the  other  hand,  if  S2  is  less  than  S-^,  turnover  can  occur  when  the  vehicle 

oscillates,  because  the  critical  slope  for  turnover  of  an  oscillated  vehicle 
will  have  already  been  surpassed  when  the  slope  is  reached  at  which 
oscillation  occurs. 

The  above  discussion  indicates  that  if  S2  is  greater  them  S^, 
turnover  will  occur  at  slope  S2,  while  if  S2  is  less  than  S-^,  turnover 
can  occur  at  slope  Sj_.  In  other  words,  the  actual  critical,  slope  for 
turnover  will  be  the  maximum  of  Si  and  S2. 

Next,  a  similar  procedure  is  used  for  oscillation  and  turnover  to 
the  right. 

Finally,  turnover  over  the  front  axle  is  obtained  in  a  manner 
similar  to  the  non-oscillating  case. 

The  equations  for  transformation  of  coordinates,  used  in  obtaining 
the  combined  CXI  for  an  oscillated  vehicle,  are  discussed  in  Appendixes  I, 

II,  and  III. 

b.  Case  II  Vehicles 


This  part  of  the  program  concerns  vehicles  with  midrange 
oscillation.  In  case  the  vehicle  has  midrange  oscillation  and  arti¬ 
culated  steering,  it  is  assumed  that  the  oscillation  joint  QJT  is 
immediately  in  front  of,  or  immediately  to  the  rear  of,  the  articulation 
joint  AJT.  For  a  vehicle  with  midrange  oscillation  and  rear  wagon  wheel 
or  Ackermann  steering,  the  oscillation  joint  OJT  can  be  anywhere  between 
the  rear  axle  and  the  front  axle. 

Once  the  computer  learns  that  the  vehicle  under  consideration 
has  midrange  oscillation,  it  goes  through  the  following  steps.  The 
computer  first  calculates  the  CG  and  weight  of  the  front  part  of  the 
vehicle,  the  CG  and  weight  of  the  rear  part  of  the  vehicle,  and  the 
combined  CG  and  votal  weight  of  the  vehicle. 

The  compute!  then  varies  the  downslope  direction  G 
(measured  counterclockwise  from  the  positive  x-axis)  from  0  degrees  to 
360  degrees,  one  degree  at  a  time.  For  each  downslope  direction  9, 


the  plane  on  which  the  vehicle  is  resting  is  tilted  in  the  direction  0 
until  at  least  one  of  the  four  wheel-loads  (as  calculated  by  the  method 
described  in  section  3),  becomes  zero.  Thus,  for  each  downs lope  direction 
0  we  obtain  the  critical  slope  angle  0.  For  0,  the  computer  prints  0,  0, 
tan(0)  in  percent,  and  the  four  wheel-loads  corresponding  to  the  downslope 
direction  0  and  the  slope  angle  0.  If,  for  a  given  0  and  0,  only  one  of 
the  upslope  wheel-loads  is  zero,  we  can  assume  that  we  have  a  situation 
where  oscillation  is  about  to  occur.  If  two  (or  three)  of  the  uphill 
wheel-loads  are  zero,  then  we  have  a  situation  where  turnover  is  about 
to  occur.  After  the  table  of  values  for  0,  0,  tan(0),  and  the  wheel - 
loads,  is  obtained,  the  table  can  be  examined  to  determine  which  down- 
slope  directions  give  oscillation  and  which  ones  give  turnover,  and  the 
corresponding  minimum  slopes  can  be  obtained  for  each  case. 

After  obtaining  the  table  of  critical  r lopes  using  the  wheel-load 
method,  the  computer  finds  critical  slopes  for  turnover  of  an  oscillated 
vehicle  by  assuming  that  the  vehicle  is  completely  oscillated  and  is 
resting  on  three  tires  with  the  fourth  tire  off  the  ground.  Four 
critical  slopes  for  turnover  are  obtained  by  assuming  that  each  of  the 
four  tires  are,  in  turn,  off  the  ground.  These  four  values  can  be 
examined  to  determine  whether  or  not  the  vehicle  will  turn  over  when 
a  slope  is  reached  such  that  oscillation  occurs. 

Example  outputs  for  both  Case  I  and  Case  II  vehicles  are  given 
Appendix  V. 

c.  General  Flow  Diagram 

The  diagram  shown  in  Figure  13  gives  a  general  idea  of  the 
main  sections  of  the  computer  program. 

d.  Input 

In  this  section  we  will  give  a  description  of  each  of  the 
21  input  cards  for  the  Static  Stability  Study  computer  program.  The 
first  card  is  a  heading  card,  and  the  second  and  third  cards  enable  the 
computer  to  decide  which  of  the  18  vehicle  data  cards  (cards  4  through 
21)  are  to  be  read. 
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Figure  13.  General  Flow  Diagram 
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(1)  Card  1  contains  the  heading,  vehicle  description,  or 
any  other  desired  information  pertaining  to  the  run.  A3 1  80  columns 
of  this  card  will  be  printed  as  the  first  line  of  output. 

(2)  Card  2  contains  (right  justified  in  columns  1  and  2) 
the  integer  NCD  -  the  number  of  vehicle  data  cards  to  be  read.  If  NCD  = 

18  (the  maximum  number  of  vehicle  data  cards),  then  card  3  must  be  omitted 
and  all  18  vehicle  data  cards  will  be  read  in  order. 

If  NCD  =  0  (i.e.,  if  card  2  is  blank),  then  the  program  will 

stop. 


(3)  If  0<NCD<18,  then  card  3  must  contain  the  integers 

n  ,  n  ,....,  nvl„_  where  n  is  right  justified  in  columns  1  and  2, 

1  ^  NCD;  1 

ri2  is  right  justified  in  columns  3  and  4,  and  iig  is  right  justified  in 

columns  5  and  6.  The  computer  will  then  read  the  vehicle  data  cards 

with  identification  numbers  n. ,  n„,...,  nVI__,  respectively. 

1  ^  NCD 


If  NCD  =  18,  card  3  must  be  omitted. 

The  18  vehicle  data  cards  contain  information  as  follows, 
parentheses  indicate  the  dimension  of  a  variable;  for  example, 
indicates  that  the  card  is  to  contain  the  three  coordinates  of 


point  FL. 


(4) 

Data 

card 

#1 

(5) 

Data 

card 

92 

(6) 

Data 

card 

93 

(7) 

Data 

card 

ft  4 

(8) 

Data 

card 

#5 

(9) 

Data 

card 

96 

(10) 

Data 

card 

ft" 

(11) 

Data 

ca  /d 

98 

(12) 

Data 

card 

.99 

(13) 

Data 

card 

910 

(14) 

Data 

card 

911 

(15) 

Data 

card 

912 

(16) 

Data 

card 

91? 

(17) 

Data 

card 

914 

(18) 

Data 

ca.rd 

915 

(19) 

Data 

card 

916 

(20) 

Data 

card 

917 

(21) 

Data 

card 

918 

The 
FL(  3) 
the 


FLO) 

FRO) 

RR(  3  ) 

RL(  3  ) 

AJT(3) 

0JT( 3 ) 

FAXO) 

CFO),  WF 
CLD(  3)  ,  WLD 
CXFO),  WXF 
CR( 3  ) ,  WR 
CRAX(3) ,  WRAX 
CXRO),  WXR 
ALPHA,  in  degrees 
OSC,  in  degrees 
SHIFT 
N STEER 
NOSC 
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On  data  cards  #1  through  # 16,  ten  columns  are  allowed  for  each 
value.  Each  value  should  contain  a  decimal  point  (on  data  cards  #1 
through  #16)  and  should  be  right  justified  in  its  respective  ten  columns. 
Hence,  30  columns  will  be  used  on  data  cards  #1  through  #7,  40  columns 
will  be  used  on  data  cards  #8  through  #13,  and  10  columns  will  be  used 
on  data  cards  #14  through  #16. 

On  data  cards  #17  and  #18,  the  integers  NSTEER  and  NOSC  are 
punched  in  column  1 . 


For  a  given  run,  the  number  of  vehicle  data  cards  present  is 
given  by  NCD  on  the  second  input  card.  The  vehicle  data  cards  whose 
ID  numbers  are  not  included  in  the  list  n.^,  n2,  .  .  • 
should  be  omitted. 


nNCD,(if  NCD  < 18) 


For  further  details  on  the  input  formats,  the  program  listing 
should  be  examined  (Appendix  IV). 

5.  Conclusions 


The  methods  discussed  in  this  report  should  be  adequate  for  study¬ 
ing  the  static  stability  properties  of  many  types  of  rough  terrain 
vehicles.  It  should  be  realized,  however,  that  the  methods  used  herein 
are  based  on  equations  for  rigid  mechanical  systems  and  cannot  be 
expected  to  give  results  that  agree  completely  since  most  vehicles  to 
be  considered  will  be  of  the  low-pressure,  rubber-tired  type.  It 
should  be  remembered  that  in  most  cases  the  center  of  a  wheel  load 
will  be  several  inches  uphill  from  the  center  of  the  corresponding 
tire  footprint.  Lack  of  knowledge  about  the  exact  location  of  the 
center  of  the  wheel  load  could  quite  easily  lead  to  computer  results. 
However,  it  is  expected  that,  through  careful  calculations  of  the 
centers  of  the  wheel  loads,  results  will  be  obtained  which  agree  quite 
closely  with  experimental  data. 

Since  the  computer  program  is  much  faster  when  the  action  line 
method  is  used  as  compared  with  the  wheel  load  method  it  was  hoped 
that  every  type  of  vehicle  could  be  analyzed  by  the  action  line  method. 
However,  it  was  discovered  that  this  method  would  not  work  for  a 
vehicle  with  midrange  oscillation;  that  is,  one  cannot  determine  the 
critical  slope  for  oscillation  to  occur  for  a  vehicle  of  this  type* 
therefore,  it  was  necessary  to  go  to  the  wheel  load  method  for  vehicles 
with  midrange  oscillation  to  determine  when  oscillation  would  occur. 
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The  present  computer  program,  which  utilizes  both  methods,  requires  about 
8,000  twenty-bit  words  of  memory.  Any  increase  in  the  size  of  the 
present  program  would  cause  the  program  to  overflow  the  memory  of  many 
of  the  smaller  computers. 

It  is  anticipated  that,  through  use  of  the  accompanying  computer 
program,  the  objective  of  obtaining  a  more  efficient  method  of  analyzing 
the  static  stability  properties  of  rough-terrain  vehicles  will  be  accom¬ 
plished.  In  addition  to  analyzing  existing  vehicles,  the  computer 
program  should  prove  quite  useful  in  designing  vehicles  by' providing  a 
raptd  means  of  comparing  the  static  stability  properties  of  two  or  more 
proposed  vehicles. 

If  it  is  decided  to  make  a  study  of  the  dynamic  case  in  which  driving 
forces,  acceleration  forces,  and  centrifugal  forces  must  be  taken  into 
account,  a  modification  of  the  wheel  load  method  will  probably  be  most 
successful  for  this  situation  since  the  wheel-load  equations  can  be 
adapted  without  tjo  much  trouble  tu  handle  additional  applied  forces. 

It  is  hoped  that  the  material  contained  in  this  report  will  prove 
useful  to  those  wishing  to  analyze  the  static  stability  properties  of 
rough-terrain  vehicles. 
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APPENDIX  I 


TRANSFORMATION  OF  COORDINATES 


The  matrix  equation[  i]jf or  transforming  the  coordinates  (p^,  p Pg) 
of  a  point  P  with  respect  to  the  x,y,z  axis  system  to  the  coordinates 

»  t  »  *  »  t  » 

(Pl»  P2»  P3'  of  the  same  point  with  respect  to  the  x  ,  y  ,  z  axis  system 
(Fig.  1.)  is  as  follows: 

(1)  P*  =  T(P-Q),  where 


Figure  1.  Transformation  of  Coordinates 


where  p^,  p2,  P3  are  the  coordinates  of  ?  with 


respect  to  the  old  axis  system; 


Q  “  J  q2  j  ,  where  qlt  q2,  q3  are  the  coordinates  of  the 
>V3, 

new  origin  with  respect  to  the  old  axis  system; 


p  '  p2  I  ,  where  p^,  p^,  p3  are  the  coordinates  of  P  with 


respect  to  the  new  axis  system; 


T 


where  tn,  tyz,  and  t]_3  are  the 


direction  cosines  of  the  x  axis  relative  to  the  x,  y,  and  z  axes 
respectively;  t21,  t22,  t23  are  the  direction  cosines  of  the  y* 
axis  relative  to  the  xf  y,  and  z  axes  respectively;  and  t31,  t32,  ■'••33 
are  the  direction  cosines  of  the  z’  axis  relative  to  the  x,  y,  and  z 
axes  respectively. 


Written 

in  scalar  form. 

equation  (1) 

becomes : 

. — 

pl  = 

tn<pi~qi> 

+ 

t12(p2'q2)  + 

t13(p3~q3) 

<  P2  = 

t2i(prqi> 

j. 

t22(p2-G-2)  + 

WW 

j>3  = 

t31(pl-ql) 

+ 

t32(p2_q2)  ~ 

t33('P3~q3^' 

These  equations  are  used  in  the  Static  Stability  Study  computer  program 
whenever  it  is  necessary  to  transform  coordinates  from  one  axis  system  to 
smother. 
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APPENDIX  IT 

ROTATION  OF  A  POINT  ABOUT  A  LINE 

In  recalculating  coordinates  in  the  Static  Stability  Study 
computer  program,  it  is  quite  often  necessary  to  rotate  a  point  Xq 
about  a  line.  Here,  we  describe  the  equations  used  to  find  the  point 
Xj  into  which  the  point  XQ  is  rotated. 

Consider,  as  in  Figure  1,  the  line  L  determined  by  the  point  J 
and  the  unit  vector  U.  Suppose  that  the  point  Xq  is  rotated  about  the 
line  L  by  an  angle  where  positive  rotation  is  taken  to  be  such  that 

the  right  hand  rule  will  give  a  direction  parallel  to  U. 


Figure  1.  Rotation  About  A  Line 
The  problem  is  to  find  Xlf  the  point  into  which  Xq  is  rotated. 

Let  P  be  the  projection  of  Xq  onto  the  line  L.  We  wish  to  find  P 
in  terms  of  known  quantities.  Since  J  and  P  are  on  the  line  L,  the 
vector *P  -  "j  must  be  parallel  to  L,  i.e., 
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Equations  (4)  and  (5)  are  used  in  subroutine  SROT  of  the  computer 
program  to  recalculate  coordinates  of  a  point  on  a  vehicle  after  part 
of  the  vehicle  is  rotated  relative  to  the  rest  of  the  vehicle. 


APPENDIX  III 


TRANSFORMATION  OF  COORDINATES  FOR  OSCILLATED 

VEHICLES 


In  computing  turnover  slopes  for  an  oscillated  vehicle,  the 
computer  program  assumes  that  the  vehicle  is  resting  on  three  tires. 

We  will  describe  the  method  used  to  transform  coordinates  so  that  the 
x,y  plane  contains  the  three  points  on  which  the  vehicle  is  assumed 
to  be  resting. 

The  problem  is  as  follows: 

Given  three  points  A,  B,  C  as  in  Fig.  1,  find  the  transformation 
equations  so  that 

(1)  the  new  origin  is  the  midpoint  of  the  segment  AB, 

(2)  the  new  y-axis  contains  the  segment  AB, 

(3)  the  new  z-axis  is  perpendicular  to  the  plane  ABC,  and 

(4)  the  new  x-coordinate  of  C  is  positive  (or  negative, 
as  desired). 


Figure  1.  Transformation  of  x,y  Plane  to  the  Plane 
Containing  the  Three  Tires  on  Which  the  Vehicle  is  Resting 
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Conditions  (1)  through  (4)  are  set  up  so  that  the  points  A,  B 
will  be  the  rear  tires  if  the  new  x-coordinate  of  C  is  to  be  positive; 
otherwise  the  points  A,  B  will  be  the  front  tires. 


Let  X,  Y,  Z  denote  the  direction  cosine  vectors  of  the  three  new 

****♦■  ► 

axes,  respectively,  relative  to  the  old  axes.  We  wish  to  find  X,  Y, 
and  2.  Let 


i 


B  -  7T,  and 

(C  -~A)  X  <B  ~~A)  =  (C  -“ A)  XYj  . 


Then  Y^  and  have  the  correct  directions  (for  the  case  shown  in 
Figure  1),  but  they  are  not  unit  vectors,  as  direction  cosine  vectors 
must  be.  Therefore,  we  set 

*Y  =  <l/|Yi/)Yi, 

(6)i^  =  a/\zx\)%,  and 

~X  ~  ~Y  X"Z. 


Then  X,  Y,  and  Z  are  the  required  direction  cosine  vectors  if  we  wish 
the  new  x-coordirate  of  C  to  be  positive. 


If  we  wish  the  new  x-coordinate  of  C  to  be  negative,  then  instead 
of  equations  (6)  we  must  use 

Y  =  (1/|Y1|)Y1, 

(7>{”z  =  and 

X  =  *?  X~$. 

^  ^ 

In  general,  once  we  have  obtained  Y^  and  Z,  from  (5),  we  can  get 
X,  Y,  and  Z  from 

^  =  d/lY^hY,, 

,  and 

"x  =  “y  x“£. 

where  B=  1  or  -1  depending  on  whether  we  wish  the  new  x-coordinate 
of  C  to  be  positive  or  negative. 
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4?  * 


Once  we  have  X,  Y,  and  Z  from  (8),  we  can  t ran si orm  any  point  using 
the  equations  described  in  Appendix  I,  where  X,  "y,  and  "z  are  the  first, 
second,  arm  third  rows  respectively  in  matrix  T,  and  Q  is  the  new  origin 

=  \  (A  +"E). 

The  methods  described  here  are  used  in  subroutine  OSN  of  the 
computer  program. 


n  o 


APPENDIX  IV 


Following  is  a  FORTRAN  listing  of  the  static  stability  study 
computer  program: 


C 


COMMON  storage  for  stability  study.  ME-010 


c 


c 

rj 


t 


C 


c 


n. 


c 


c 


COMMON  AJT13],  AL113],  AL2t3l.  CF  {  3  J  ,  CG  [  3  }  ,  CLD  f  3  I  »  CR  (3  J  .  Cf'AX  t  3  ! 
1.  UXMA1.  CXR  !  3  )  •  F  (  4  1  »  F  A  y  J  3  )  .  Fl(3I.  FlP[3).  Fk{3I,  FrtP(3l 
COMMON  0 JP I  3  J  »  OjTIJI,  RI9.3J.  RL  t  3  J  #  RLP(3i.  RR  t  3 1  •  RRP131,  U(3). 
l<(vl.  vp ( 9 ! ,  XPP ( 9  1  »  Y{9],  YP191.  YPP[9)>  Z  l  9  )  .  ZPI9).  ZPP[9) 

10^'IN  ALPHA,  DE'j  ,  UA,  CSC,  PHI,  SHIFT,  T,  IH,  TPH ,  U,  W1  ,  w?.  wF , 
1WL-J.  WP,  WHAX.  wXF,  WXR 

common  1REAU,  NOSC,  K'STEER 


cj MENS! UN  uF  P I  3 ! . CL  UP ( 3 ) , COSC t 3 ) . CST [ 3 1 . CTnO I  3 1 , CTO [ 3 J , CX  l  3  1  , 
3  .Xf  P(31  ,F  AXP  m  ,FC  R]  .-UX  t  3  !  ,UZ  [  31 
■'  I  MENS  I  ON  NtlM 

READ  ANU  PRIM  HEADING 
mii  4  E  4  0  i 
3  ►  OHM  AT  (BOH 

1  1 
^HiKT  h 

R  I  M  T  1 

REAL)  NCL)  =  NU  OF  DATA  CARDS 
■<b  * H  ?  .  NC i) 

'  U«MAT  ( 10  IP,  in  I  >  1 

!  F  •  A.CO  1  lijuun  .lOf'OO,  1  on 
in-.  iFiNcn-id]  3no,?o»,2nn 
'U  *  'o  ?m  i  =i.ip 
'  i'i  \  [  I  I  s  i 

■?u  to  Ann 

.Mjtl  4=AD  ?,  INI  I  1  .  1=1, NCrn 

RP  AD  AND  PRINT  VFHICLE  DATA 
**•)•»  JR  i  mT  3 

*>  *■  OMMAT  i ///.  30  *  •  l^HVtH  I  CL  It  DATA,///] 

■  nj  M.in  I  =  l , ncl) 

!  R~  a 0  =  N(I) 

•'All.  RP  AO 
•’l,  ..  MIL  PRINT 
J  R  I  'J  T  S 
*7  >  OKNAT  I  1  H  1  ] 

ir.i  =  3.141b9?6b  /  IPO. 

ALCRR  =  ALpHA  *  OtG 

RECALCULATE  COORDINATES 
0  7  1  I  =  0  . 

07 i /  ■  =  n  . 
i  1 7  t  3  1  =  1 . 

••0  MO  I  =  1,3 
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*  "  I 


n  o  o  non 


m  n 
C 

C 

c 

60  0 

c 

c 

/on 


>10 

c 

c 

600 

900 

10  0  0 
6 


FLPCI]  s  FLU) 

FRPI11  a  FR( 11 
RRP [  I  1  «  RR  t  !  1 
oi  p  r  i  l  i  oi  m 

’  w  *•#  •  »  *—  •  »  . 

l  <  o  f  oi  •  ri  roi  -  cutr r 

1  W.'  I  <  —  •  ••  »  fc  I  «-»■••»» 

F  R  P  (  2  I  *  FR123  +  SHIFT 
RRP I ? ]  *  RR [ 2  1  *  SHIFT 
: < L P  t  2 1  s  RU21  -  SHIFT 

GO  TO  [600. 700,600] .NSTEER 

CALL  SROTIAJT.UZ.ALPHR.FAX.FAXPI 
CALL  SROTIAJT.UZ.ALPHR.CF  ,CFP  J 
CALL  SROTt AJT.UZ. ALPHR.CLD,CLDP3 
CALL  SROTIAJT.UZ.ALPHR.CXF.CXFPl 

CALL  SR0T [AJT.UZ* ALPHR,f LP,FLP  ) 

CALL  SROT [AJT.UZ* AlPHR*FRP,FRP  ] 

GO  TO  800 
U  Z l 3 ]  =  -UZ 13  3 

CALL  SROT[CHAX*UZ, ALPHR.RRP.RRP] 

CALL  SROTICRAX.UZ, AlPhR,RLP#RLP1 
DO  71 0  I  =  1,3 
CFP[1)  =  CM  I  1 
C X I-  P  [  I  ]  =  CXF  [  I  ] 

CL ‘JP  i  I  1  =  CLOf  i  1 

CHECK  TYPE  OF  OSCILLATION 

GO  TO  (900.900.60001.  NOSC 
DO  5000  NDS  s  1.3 
GO  TO  (1000.1300.4000].  NDS 
PRINT  6 

FORMAT!//. 50X.24HDOWNSLOPE  IS  TO  THE  LEFT 
NON-OSCILLATED  TURNOVER  CALCULATIONS 
COMPUTfc  COMBINED  weight 


WTNO  =  WRAX  ♦  wR  ♦  WXR  ♦  WF  ♦  WXF  ♦  WLD 
COMPUTE  COMBINED  CG 


DO  1010  I  =  1,3 

CTNOtU  =  l  CRAX 1 1  I  *WR  AX  ♦  CR(U*WR  ♦  CXR[J]*WXR  ♦  CFP(I]*wF 
1  ♦  CXFP I  I ] *WXF  «■  CLDP(I3*WLD  1  /  WTNO 

mill  CGI  I  ]  =  CTNOl  1  ] 

GO  TO  1300 
1100  PRINT  7 

7  KORMATI//.DOX.25HDOWNSLOP6  'S  TO  THE  RIGHT  .//l 

no  in n  i  =  1,3 

1110  CG  t  I  ]  =  CTNOtl J 
13011  PRINT  8 

ri  FORMAT [//.30X.42HSIUESLOPE  TURNOVER  ASSUMING  NO  OSCILLATION  .//] 


C  CALCULATE  ACTION  l INE 

co  io  1 1 4 u o » i 5 o n j ,  nds 

1 A  U  i)  ')0  1410  I  =  1,3 

ALL  I  ,  )  =  RLPtll 
1*10  fits r  I  j  =  FLP r I  1 

c 

■iO  TO  1550 

c 

15JII  no  1510  I  =  1,3 
ah  [  1  ]  =  RHP  t  I  ] 
lult)  4L<-  {  I  1  s  FPPt  I  1 

C 

1TJ5.I  CALL  SLOPE 

CO  TO  15000,16001 ,  NOSC 
C 

LUJ.I  PRINT  0 

7  FijRMAT[//,30x.58HSIOeSLOPE  AT  WHICH  COMPLETE  OSCILLATION  IS  LIKELV 
1  TO  OCCUR  ,//] 

C 

r  COMMUTE  CG  and  ACTION-LINE  for  vehicle  minus  rear  axle 

C 

WQSC.  =  WTNO  -  WRAX 
00  161 n  I  =  1,3 

c 

coscm  =  i  criiuwh  +  c<r:ii*wxr  ♦  cfpuj*wf  *  cxfpiii*wxf 

1  ♦  CLOP  1  I  1 *WLO  1  /  wnsc 

o 

CG -11  =  COSCl  I  1 

1*  1  !  a L  1 1  1  !  S  CKAXI  I  1 

c 

~0  TO  13700,1800),  NDS 

c 

i/o»  co  i/in  i  =  1,3 

1  / 1 »»  ALc  III  =  FLPt  I  ] 

r 

GO  TO  1 P50 
C 

LMin  > )  0  1610  I  =  1,3 
mlO  AL'l  I  1  =  fhp  r  I  ] 

C 

i-ou  caul  slope 

C  CALCjtATIONS  for  oscillated  TURNOVER 

PRINT  10 

1!>  -06MAT[//,3UX,51HTURNCVER  ASSUMING  COMPLETE  OSCILLATION  HAS  OCCURH 
1  Ml  .//] 

C 

TA'J  =  r  t-1.  ]**{NDS*1  1  3  *  OSC  *  DEG 
C 

ox  .  1  1  =  -1  • 

"  X  !  ?  1  =  >1  . 

1 1 X  t  3 1  =  n. 

C 
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CALL  SROT [LRAX.UX.TAU.COSC,  CST] 

C 

GQ  T 0  (1900,9000  ).  NL1S 

C 

1900  CALL  SROT(CRAX,UX,TAU,FL,FC) 

GO  TO  2010 
C 

2 o 0 0  CALL  SROTtCRAX.Ux.TAU.FR.FC) 

C 

2010  MEL  =  1. 

C 

CALL  OSM  RR » RL » F  C  .  DEI  .CRAX.CX) 

CAul  OSNCRR,RL,FC,DEL.CST,CST) 

CALL  OSN(RR,RL,FC.rjEL.FC,  Al2  ) 

C 

00  2020  I  =  1.3 

CTO(I)  =  (  CSTIIUROSC  ♦  CXll)*WRAX  )  /  WTNQ 
21*20  CGi  H  =  CTO  t  I  ) 

C 

GO  TO  (2100,2200),  NOS 
C 

>'1  00  1)0  2110  I  =  1.3 
2110  AL1  (  I  )  =  RL  (  I  1 
C 

GO  TO  3000 
C 

220*1  nn  2210  I  =i,3 
<  <r  1  '■  Al.  tin  =  RKl!) 

C 

3 (■')•*•  CALl  SLOPE 
C 

GO  TO  S000 

CALCULATIONS  for  turnover  over  front  AXLE 
«O0r*  PRiNT  11 

11  S-'ORMAT  (//.t)OX,2^HOOwNSLOi3fc  IS  TO  THE  FRONT  ,////,  30X ,  2BHTURN0VEP  0 
lvE«  THE  FRONT  AXLE  ,//  ) 

C 

C  CALCULATE  ACTION  LINE  and  combined  cg 

c 

GO  TO  [4100.4200,4200!.  NSTEER 
C 

MOO  CALL  SR0T(AJT,UZ»AL?HR,FL.  AL1) 

Call  SPOT  [  A  JT,  U?..  ALPHR.FR,  AL  2) 

C 

GO  TO  4300 
C 

4  2  1)  *■  uo  4210  I  =  1.3 
ALIM!  =  FLII  1 
42  '  )  ALUM  =  FR  [  I  1 
r 

•*3  0  0  00  4310  I  =1,3 
4310  Coll)  =  CTNO ( I  ) 


o  o  o  r> 


C 

CALL  SLOPE 
C 

C  T HE  NEXT  STATEMENT  ? S  THE  END  OP  THE  NDS  LOOP 

500!)  continue 

C 

CO  TO  9000 

CALCULATIONS  fop  a  vehicle  WITH  midrange  oscillation,  using 

THE  WHEEL-LOAD  METHOD 

e-jr  I  IKINSTEER  -  1!  6300,61  00,6300 
Oluu  I F i OJT  f  2  J  -  AJTI2)  I  6300,6200,6200 
C  CALCULATE  UNIT  VECTOR  PARALLEL  TO  OSCILLATION  AXIS 

e2U0  DO  6210  I  =  1,3 

62iii  uin  s  faxpi 1 1  -  ajti n 

c 

Au  =  SORTF  (  U ( 1 1 **2  ♦  U I  2 1 **2  ♦  UE31**2  I 
C 

DO  6220  I  =1,3 
U  [  i  1  =  lit  I  I  /  AU 
C 

CAuL  SPOT [AJT.UZ, ALPhr.OJT.OJP) 

C 

)  T U  64  0  0 

r* 

U 

6  .i  it  -i  U  !  L  )  -  1  . 

•  i[^i  =  n . 

IJ'AI  =  0. 

c 

■iQ  ftUH  I  =  1,3 
nil  I!  0  JM  J  =  OJTm 

£ 

fc h  >•  *1  =  hF  *  WLQ  ♦  **XF 
•S?  =  WPAX  ♦  HP  ♦  WXH 
W  =  W 1  ♦  W  2 

C 

C  HAKE  CALCULATIONS  PRELIMINARY  TO  CALL  TO  WHEEL-LOAD  SUBROUTINE 

C 

DO  70  on  J  =  1,3 
C 

R  f  1  , J ]  =  FLPIJ] 

R  l  k  » .1 J  =  FRPI  J) 

r  [  i ,  j  i  =  rhp  r  j ) 

PS-*,JI  =  RLP  l  J  ] 

=  (  CFP  [  J 1  *wF  ♦  ClDPIJJ*WLD  ♦  CXFPlj]*wXF  ]  /  W1 
HIo.J]  =  (  CRAXI Ji*wRAX  *  CR  E  J I *WR  ♦  CXREJ]*WXR  ]  /  W2 
Rf/.J]  =  (  R 1 5 , J 1 *W1  ♦  H  E  6 , J I *W2  i  /  W 
H i i , J J  =  O jP  E  J I 
RfY.Jl  =  UIJI 
•  ti 0 "  CONTINUE 

C  PND  OF  J  LOOP 

uo  ;-ino  I  =  1,9 
X  f  11  =  R  f I » 1 1 


U  0 


ooooo  oooo  a  orto 


front  AND  REAR  WEIGHTS, 


PRINT  12,  X'51.Yt5],2t5I,Wl#Xl6J*YC6j,2l6J,W2,X(7J#Yf7l,7l7J,W 

12  FORMAT  [ 1  OX » 1H [ ,F9.1»1H» ,F9.1,1H, »F9«1,1HJ  »11H  *  FRONT  CG,F20.1,l5H 
1  s  FRONT  WEIGHT  , // , 1  OX , 1H ( , F9 . 1 , 1H. , F9 . 1 , 1H , , f9 . 1 , 1H ) , 1 3H  *  REAR 
2CG  .F2l.-1.14M  =  REAR  WEIGHT  , //, 10X, 1H ( , F9 . 1 , 1H, , F9 . 1 , 1 H, , F9 . 1 , 1H I 
3,14H  =  COMBINED  CG  .F17,1,15H  *  TOTAL  WEIGHT  .//  J 

l 

MAKE  WHEEL-LOAD  CALCULATIONS.  VARY  ORIENTATIONS,  AND  PRINT 
CORRESPONDING  CRITICAL  SLOPES  AND  WHEEL  LOADS. 


CAUL  MIDOSC 

FIND  CRITICAL  SLOPES  AND  DOWNSLOPE  DIRECTIONS  FOR  OSCILLATED 
TURNOVER,  FOR  THE  VEHICLE  RESTING  ON  THREE  TIRES,  WITH  THE  4TH 
TIRE  OFF  THE  GROUND.  (FOUR  CASES! 


CALL  MDOTVR 
C 

V 0  0 0  GO  TO  50 
10000  PRINT  13 

li  FORMAT (/////, 50X,l5HTHATS  ALL  FOLKS  J 

C 

call  exit 

END 


u 
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READ 


this  subroutine  reads  the  data  card  with  id  number  iread.  and 

THEN  RETURNS  TO  THE  MAIN  PROGRAM . 

SUHROUTINE  read 

COMMON  STORAGE  FOR  STABILITY  STUDY,  ME-010 

COMMON  A JT ( 3 ] ,  AL1 1  31 »  AL2I31,  CF131,  CG131*  CLDt3)»  CP ( 3 1 , CRAX t 3 ] 
1,  CXFI31,  CXR [ 3 1 ,  F  [  4  1 ,  FAX131,  FLI31,  FlP(3),  FRI3),  FRP13) 

COMMON  0 JP  131,  0JT131.  R ( 9 , 3 ) ,  RL131,  RLPI31,  RRI31,  RPP(3)»  U  t  3 ) , 
lx[v]»  XP I  9 i ,  XPP  t  9  1  ,  Y  t  9 1 ,  YP191,  YPP  t  9 1 ,  Z  t  9 1 »  7P  t  9 1  ,  ZPPtVl 

COMMON  ALPHA,  DEG,  OA,  OSC,  PHI,  SHIFT,  T,  Th,  TPH,  W,  Wl»  w2,  WF . 
lWL 0 ,  WR *  WRAX.  WXF.  HXR 

COMMON  I  RE AO ,  NOSC,  NSTEER 


1  >1  0 
2  0  0 


1 

> 

;s 

4 

*5 

/ 


l'J 

11 

12 


15 


1  a 

l'J 
16 
1  I 

n 


FORMAT (flFlO.O  1 
FORMAT (Ill 

IREAD  =  THE  ID  NUMBER  OF  THE  DATA  CARD  WE  WANT  TO  READ. 
IR  =  IREAD 

CO  TO  [  1,2, 3, 4, 5, 6, /, 8, 9, 10*11, 12,13, T4, 15, 16, 17. 18  1  »  I* 
READ  100,  fl 
RE l URN 

READ  inn,  FR 
RE  I  URN 

REAP  100.  RR 

RE i URN 

REAP  100.  RL 
re  Turn 

read  inn,  a jt 

RE  I  URN 

READ  100,  0 JT 
Rt • URN 

READ  100.  FAX 
uF I  URN 

RFAp  100,  CF,  WF 
RE  I  URN 

REAP  100.  CLD  ,WLD 
RE  » URN 

RF  AO  100,  CXF ,  WXF 
RE i URN 

RE  AO  100.  CR,  W« 

RE  I  URN 

READ  inn,  URAX,  MR  AX 
RE  I  URN 

REaD  100.  CXR,  WXR 
RE  !  URN 

READ  IOC,  ALPHA 
RE i URN 

read  ion,  OSC 

RE  TURN 

READ  100.  SHIFT 
URN 

READ  200,  NSTEER 

q c  ]  Mt|w 

hFAD  200.  NOSC 
RE i URN 
c-NO 
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PRINT 


THIS  SUBROUTINE  PRINTS  THE  DATA  CARD  WITH  ID  NUMBER  {READ,  AND 
THEN  RETURNS  TO  TmE  MAIN  PROGRAM, 

1 

SUBROUTINE  PRINT 

COMMON  STORAGE  FOR  STABILITY  STUDY.  ME-010 

COMMON  A JT 1 3  j  *  All 1 3 ] ,  AL2I3J,  CFI3).  CG  t  3 1 »  CLD ( 31 .  CR t 3 J , CRAX 1 3 1 
1.  CXFI3J.  CXR ( 3 ) .  F  1 4 ] .  FAXI3I,  FL  I  3 1 »  FLP[3I.  FRI3J.  FRP(3I 
COMMON  0  JP 13  3*  0JTI31,  RI9.3I.  RU3J.  RLPI3I,  RRI3I.  RRP 13).  U 1 3 1  . 
lxm.  xP  1 9 1 »  XPP  I  9 1  *  Y  t  9 1 .  VP  1 91  •  YPP  191.  Z19I,  ZPI9).  ZFP  t  9  I 

COMMON  ALPHA,  DEG,  OA,  OSC.  PHI,  SHIFT,  T,  TH»  TPH,  W»  Hl>  W2»  W?  , 
1WL  L) »  HR,  WRAX,  WXF »  WXR 

COMMON  IRF.AD,  NOSC,  NSTEER 
C 
C 

IR  s  IREAD 

GO  TO  I  1.2,3,4,5,6,7,8,9.10,11,12,13,14,15,16,17.18  1  ,  IR 

C 

1  PRINT  100,  FL 
RETURN 

2  PRINT  ?00,  FR 

return 

3  PRINT  300,  RR 
R61 URN 

4  PRINT  400,  RL 
RE  TURN 

5  PRINT  500.  A JT 

return 

6  PRINT  600.  OJT 
RF I  URN 

7  print  700,  fax 

RETURN 

5  PRINT  800,  C i’»  WF 
RETURN 

9  PRINT  900,  ClD,  WlO 
RE  fURN 

10  PRINT  1000.  CXF.  WXF 
RE  I  URN 

11  PRInT  1100.  CR.  HR 

return 

12  PRINT  1200*  CRAX.  WRAX 

return 

13  PRINT  1300.  CXR.  WXR 
RETURN 

14  PRINT  1400.  ALPHA 

return 

15  PRINT  1500.  OSC 
RETURN 

16  PRINT  1600.  SHIFT 


RE I URN 

1/  PRINT  1700*  NSTEER 
RE  t  URN 

1H  PRINT  1MOO.  NOSC 
RETURN 

lOH  FORHATtl0X»lMt,F10.1.1H*,Fl0.1.1H,,Fl0.1»lH)#5X,3H*FL  1 
?UU  FO«MAT(10X»1H[ »F10.1»1H» *F\0.1*1H* *F10.1»1HJ »5X,3H*FR  1 
400  FONHAT[10X*lHl»F10.1tlH».Fi0.1»lH,.Fin.l*lHj.5X,3H*RR  1 
400  FO«MATtl0X.lH[,Fl0.1.1H..Fl0.1.1H,,Fl0.1»lHJ#5X,3H*RL  1 
t>0li  FORMAT  110X»1H(  ,F1Q  .I.1H.  ,FiO  .1.1H*  »F10.1»1M1  »5X,4H*AJT1 
ftOO  FORMAT [10X»1H( .F10.1.1H. »F10.1»1H# #F10.1»1H) »5X,4HsOJT l 
/On  FORMATltGX.lHt  »F10.1»1H, »Fl0.1.1H. »F10.1»1H1 #5X»4H*FAX1 
ttliO  FORMATflOX*lH(,F10.1#lH# .FlO.l.lW. »F10.1»1HJ .5X,3H*CF,F2?.1»4H  «WF 
1  1 

9  0  0  F  ORMAT  f 1 0 X » 1H I . Fl 0 . 1 » 1 H * . F 1 0 . 1 , 1H *  , Fl 0 . 1 » 1H ] . 5X , 4HsCLD , F 21 . 1  * 5H  «W 
ltDl 

1  000  FORMAT  11  OX. lHt  »F10 . l.JH,,Fl0.1,lH,,Fl0.1#lH)»5X,4HsCXF.F21.1,5H  *W 
1 XF  * 

1100  FORMATt 10X.1HI .F10.1.1H, .FlO.l.lH, ,F10.1»1H1 »5X.3HsCR»F22.1»4h  *WR 
1 J 

12  01.  FOHkAT  IIOX.IH  l  »F10 .1.1W*  »FlO  .F10.1.1H)  *5X,5H»CRAX»F20 .1 ,6H  s 

1WRAX  ] 

13  0  0  FORMAT llOX*lH[.F10.i.lH,,Fl0.1»lM.»Fl0.1.1H] » 5X , 4H*CXR , F?1 . 1 , 5H  *W 

1  X  R  i 

14(10  FORMATIF20.2.35H  =AUPRA  { DEG  1  s  ANGLE  OF  LEFT  STEER  I 
ibOO  F ORMAT [F20 -2.32H  =OSC  f  DEG] =MAX  OSCILLATION  ANGLE  ] 

1600  FORMATIF20.2.39H  s  SHIFT,  MEASURED  FROM  CENTER  OF  TREAD  ] 

1 / 0 u  FORMAT (110* b2H  sSTEER  TYPE [ 1* ART ! CULATED, 2  =  REaR  WAGON, 3=ACKERmANN] 
1  1 

leOU  FORMAT  I I10.59H  =OSClLLATION  TYPE  1 1  =  NQ  0SC»2=REAR  AXLE  0SC,3  =  M1DRAN 
1GE  0SC1  1 

ENU 
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SUHPOUT I NE  SROT(XJ,U,OGA,X0.X13 

CONS  I OtR  THg  LINE  L  DETERMINED  BY  THE  POINT  XJ  AND  THE  UNIT 
VECTOR  U.  SUPPOSE  THE  POINT  XO  IS  ROTATED  ABOUT  THE  LINE  L  ThRu 
AN  ANGLE  OGA,  WHERE  POSITIVE  ROTATION  IS  TAKEN  TO  BE  SUCH  THAT 
THt  RIGHT-HAND  RULE  WILL  GIVE  A  DIRECTION  PARALLEL  TO  U.  ThIS 
SUBROUTINE  FINDS  XI.  THE  POINT  INTO  WHICH  XO  IS  ROTATED. 
SUBROUTINE  SROT [Xj.U.OGA.XO.XlI 
DIMENSION  CPl33*P(3}»Ut3]»XQ!3l#Xlt31#XJf3] 

t 

ANALYSIS 

AK  =  [XO-XJI.U  s  SCALAR  PRODUCT 
P  =  XJ  ♦  IAK]*U  =  VECTOR 

XI  =  COS  t  OGA  J  * ( X  0*P 1  *  S 1 N I OGA ] t  U  X  [XO-PJ)  ♦  P  *  VECTOR 


AK  =  0. 

COMPUTE  AK 
00  100  I  s  1.3 

ion  ak  =  ak  ♦  txnm  -  xjrm*utn 
compute  p 
DO  POO  I  =  1.3 
poo  pm  =  x j 1 1 3  ♦  ax #  um 

COMPUTE  XI 

CP  =  CROSS  PRODUCT  OF  u  AND  XO-P 

CPU)  =  Ut«cfJ*lX0t31-P(3J  ]  -  Ul3J*IX0t2J-Pt2)J 
CPi?J  =  U  1 3 1  •  [  X  0  HJ-PilH  -  U(lI*lX0t3J-Pt3n 
c p 1 3 s  uti)*[xo(2i-P{?n  -  ui2i*txofij-Ptm 

COMPUTE  final  RESULT 
DO  3'J  0  I  s  1.3 

300  XltU  s  COSF[QGA!*tXOtn-Pm  J  ♦  S  INF  [  ”>GA  I  *CP  l  I  3  ♦  rill 
C 

RETjRN 

Emu 


C  SUdROUT 1  NS  05NtA»8»C,D!rL»P»PP) 

C 

C  GIVEN  A  TRIANGLE  ARC  /  0  A  POINT  P,  WE  WISH  TO  F I ND  THE 

C  COORDINATES  OF  THE  POINT  P  RELATIVE  TO  A  NEW  AXIS  SYSTEM  HAVING 

C  THc  PROPERTIES 

C  (II  THE  NEW  ORIGIN  IS  THE  MIDPOINT  OF  SEGMENT  AS 

c  [21  the  new  y-axis  has  the  direction  of  the  vector  ab 

C  131  THE  NEW  2-AXIS  IS  PERPENDICULAR  TO  PLANE  ABC 

C  14]  THE  NEW  X-COORDINATE  OF  C  HAS  THE  SAME  SIGN  AS  DEL ,  WHERE 

C  DEL  EQUALS  M  OR  -1. 

C 

C  ANALYSIS.  FIRST  FIND  THE  DIRECTION  COSINE 

C  VECTORS  OF  The  NEW  AXES  RELATIVE  TO  THE  OLD.  THEN  TRANSFORM  P, 

C  USING  THE  USUAL  ORTHOGONAL  MATRIX  OBTAINED  FROM  THE  DIRECTION 

C  COSINE  VECTORS. 

L 

SUHROUTINE  QSNfA.B.C.DEL.P.PP) 

DIMENSION  A[?l,B(3)«C[3)»Pl31»PPt3]*PQl3)#Xl3]#Yl31»Yl(3]»Zl3]» 

1  21131 

C  FIND  THE  3  DIRECTION  COSINE  VECTORS 

C 

no  inn  i  =  1,3 
1U1I  Y3  i  n  =  sin  -  All) 
c 

C  71  =  E  C- A  1  CROSS  Y 1 

71U1  =  t  C 1 2 1 -A  1 2 1  1 *Y1 1  3 1  -  ICI31-AI3J  )  *Y1 [ 2 1 
71-21  =  IC»3]-A13) 1*Yltl)  -  tC(ll-Atl) ]*Y113] 

7H31  =  tClll-All  J  1  *Y1 1 2 1  -  l  C  1 2 1  -  A  t  2  ]  ]  *  Yl  r  1  ! 

C 


C  NORMALIZE  Yl  AND  Zl»  AND  CHOOSE  CORRECT  DIRECTION  FOR  7-AXIS. 


ay  =  SQPTFl  Yl ! 1 1 **2  *  Yl 1 21 **2  *  Y1I3)**2  1 
A7  s  SOPTFl  Zlll)**2  ♦  Z1 1 21 **  2  *  Z1I3)**2  1 
C 

DO  200  1  =  1.3 
Y 1 i ]  =  Yl II)  /  AY 
.'OH  7  f  1  1  =  GEL*Zl  f  I  1  /  A7 
C 

C  CALCULATE  X  «  Y  CROSS  7  - 

C 


V[  11  =  Y 121*7 131  -  Y ( 3 1 *Z I  2 1 
X  E  ^  1  =  Y[3)*Zni  *  YtlJ*Zi31 
xm  =  y 1 1  j *2 r 2 1  -  y 1 2 1 *z 1 1 ) 
c 

C  TRANSFORM  P  TO  GFT  THE  COORDINATES  OF  PP. 

C 

DO  300  I  S  1.3 

3 u  0  °o  1 1 )  =  n . 

r 


00  400  I  =  1,3 

r  =  pm  -  .5*1 

PQil 1  S  POll]  ♦ 
PQi?J  =  P0I2]  ♦ 
a 0 0  P Q i 3  1  =  PQ13]  ♦ 


MD^BIII 
x  m*R 
y  m  *r 

zm*R 


i 


00  500  I  =  i»3 
SOU  DP  t I  1  =  PQ  II  J 
«E  turn 
fcNU 
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subroutine  slope 

GIVEN  THE  END-POINTS  OF  AN  ACTION  LINE,  AL1  AND  AL?,  AND  A 
CENTER  OF  GRAVITY,  CG.  THIS  SUBROUTINE  CALCULATES  THE  CRITICAL 

Slope  phi,  the  percent  of  slope  tph  [equals  igo*t animhi n »  and  the 

CORRESPONDING  DIRECTION  OF  DOWNSLOPE  OA  [MEASURED  FROM  THE 
K>5iTivE  x-Axisi. 

subroutine  slope 

COMMON  STORAGE  FOR  STABILITY  STUDY,  KE-010 

COMMON  A JT l 3 ) ,  AL1 1 3 ) ,  AL2I33,  CF133.  CG[3],  CLD[33,  CP t 3 ) . CR AX t 3 3 
1,  LXFI33,  CXR I  3  3  »  F[4),  FAX  131,  FL  t  3 1 ,  F  LP  t  3  I ,  FR[3),  FRPt3J 
COMMON  0JPI31,  0 JT { 3 ) ,  R  [  9 , 3  J ,  RLC33,  RLP[31.  RR133,  RRP[3J,  U[33, 
IX ( y  3 ,  XP l 9  * ,  XPP I  9  3 ,  Y  t  9 1 ,  YPI9J,  YPP [ 9  J ,  Z l 9  3 ,  2PC93,  ZPPt9J 

COMMON  ALPHA,  DEG,  OA,  OSC,  PHI,  SHIFT,  T.  TH,  TPH,  W*  HI.  W2»  WF, 
1 W L U ,  HR,  WRAX,  WXF  *  WXR 

COMMON  IREAD,  NOSC,  NSTEER 

D I  MENS  I  ON  DA [ 3  J , DCA [ 3  3 
IF  INOSC-33  100.200.100 
ion  PRINT  1.  CG.AL1.AL2 

1  FOWMAT[/.lH{  ,F9.1,1H,,F9.1,1’4,,F9.1,1H)  ,3H*CG»6X,1H[,F9.1»1H,  .F9.1 
1,1H, »F9.1»1H],4H«AL1,5X,1H[,F9.1,1H,,F9.1,1H,,F9.1,1HJ ,4H=AL2.//  3 

200  PRINT  2 

2  F0HMATi20X»l4HCRIT!CAl  SLOPE . 16X . 42HD0HNSL0PE  D I RECT I  ON [ ANGLE  FROM 
1  PUS  X-AXIS]  J 

CALCULATE  T 
UO  210  I=l»3 
OC  A  [  I  ]  s  CGin  -  AL1  (  I  3 
211)  DA  i  I  1  =  AL2  [  I  3  -  AL1 1  I  3 

A  s  DC A [ 3  3 
B  =  -D A [ 3  3 

C  =  DC A [ 1 3  **2  ♦  DC A ( 2 1 **2  ♦  DCA(33**2 

D  =  (-2.3*1  DA ( 1 3 *DC A ( 1 3  ♦  DA [ 2 1 *DCA [ 2 3  ♦  DA [ 3 3 *DCA i 3 3  3 
F  =  DA [ 1 3  **2  ♦  DA [ 2 ) **2  *  DA[3J**2 

T  =  I A*D  -  2 . *B*C  3  /  I B*D  -  2.#A*E3 

CPH  =  (A  ♦  B*T  3  /  SQRTFIC  ♦  D*T  ♦  E*T*T) 

TPM  r  SQRTFl  ABSFtl.  -  CPH*CPH]  J  /  CPH 

CALCULATE  PHI 
PHI  =  ATANf  [TPH  3  /  DEG 

calculate  percent  SLOPE 

TPH  =  TPH  *  100. 

CHECK  TO  SEE  IF  CG  IS  DIRECTLY  OVER  ACTION  LINE 
IF!  TPH  -  l.E-5  3  300,400,400 
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30U  PRINT  3#  PHI 

3  FORMATl/,F30.2,2lX,33HNONE-CG  DIRECTLY  OVER  ACTION  LINF.////1 
RE  1  URN 

400  XT  =  AlllU  ♦  T*DaI11  -  CGI1J 
yT  =  AL i 1 2  3  -  T-OAJ25  -  CG[2J 

CALCULATE  DOWNSLOPE  DIRECTION 
4  OUAD  ARCTAN  ROUTINE 
IF  i  XT  ]  430.420,410 
410  QA  =  0. 

CO  TO  440 

42'i  OA  =  S  I GNF  t  90  .  ,  YT  ] 

CO  TO  450 

430  OA  =  S I 6NF  1 180 . , YT 1 
"4!!  OA  =  OA  1  ATANF I  YT/XT 1  /  DEG 
450  CONTINUE 
C  END  OF  4  QUAD  ARCTAN 

C  SEE  IF  THE  VEHICLE  HAS  MIDRANGE  OSCILLATION 

IF INOSC  -  33  500,600,500 
5U'i  PRINT  4,  PHI,  TPH,  OA,  T 

4  FORMAT I/,F20.?,6H  DEG  «»F9.2,4H  PCT.F21.2.4H  DEG , F51 . 5 , 2H  =  T , //// } 

G 

600  PE  I  URN 
t-ND 


ooooo  oooooooo 


SUBROUTINE  !  DOSC 


THIS  SUBROUTINE  VARIES  THE  DOWNSLOPE  DIRECTION  THROUGH  360 

ncfiricce  Aiuit  rno  c  *  ru  nftuuot  adc  imoC^t  mu  rMCiStiocn  rntiK'TFRGi  nCK« 

UU  «nu  u  j  »  w  r  \jr\  l  w«'i  3WW  iz.  1/  t  ■>  ww  i  »  t  •>  -  -  ■*  —  —  ** 

WISE  FROM  THE  POSITIVE  X-AXIS)  TILTS  THE  GROUND  PLANE  IN  THAT 


r»  t  or^r  »  ftw 
U  l  r%  r  O  !  iun 


Illicit 
1  1L 


OAtini  Ttfti.  A r 

vurvui  »  Awrv  w  r 


tkiPTABTi  1TV 
i  '1  ^  i  *•  o  i  \>  •  *  • 


ocirucn . 
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WHEEL -L 0 ADS  ARE  PRINTED  FOR  EACH  CRITICAL  SLOPE, 


SUBROUTINE  M I DOSC 


COMMON  STORAGE  FOR  STABILITY  STUDY,  HE-010 

COMMON  A JT 1 3 ] ,  AL1I3I,  AL2I3J,  CFI3I,  CGI3),  CLD  t  3 1 ,  CR ( 3 1 . CRAX I  3 3 
1,  C  x  F [ 3  3 ,  CXRI31.  F  1 4 ] •  F  AX ( 3  3  *  FLI3],  FLPI31,  FR(3J,  FRPI3) 

COMMON  0 JP 131,  OJTI33,  RI9.33.  RLI33,  RLP I  3  3  »  RRI3),  RRPt33,  U  t  3  3 . 
1X1*3.  XPI93,  XPPI93,  Y I  9  3 ,  YP ( 9  3 ,  YPPI93,  Zt91,  7Pt93,  ZPPI91 
C 

COMMON  ALPH,',  DEG.  OA,  OSC.  PHI,  SHIFT.  T.  TH,  TPH ,  W,  Ml,  W2.  WF , 
1WI.O,  WR,  WRAX,  WXF »  WXR 

C 

COMMON  IREAD,  nosc,  nsteer 

c 

PRINT  SO 

5D  FORMAT!///. 10X.103HFOLLOWING  IS  A  TABLE  GIVING  SLOPES  aT  WHICH  INS 
1TAHILITY  OCCURS  FOR  A  VEHICLE  WITH  MIDRANGE  OSC I LL AT  I  ON . , / . 10 X » 
271HTHETA  =  DOWNSLOPE  DIRECTION  I  MEASURED  COUNTERCLOCKWISE  FROM  POS 
3  X-AXIS1  ./.10X.77HPHI  ■  ANGLE  THRU  WHICH  GROUND  PLANE  MUST  BE  TI 

4LTED  FOR  INSTABILITY  TO  OCCUR »/»  10* .  27HSLOPE  =  T  AN  [  PH  I  ]  IN  PERCENT 
5.1UX.56HFCI3  =  WHEEL  LOAD  AT  TILT  ANGLE  PHI.  AS  INDICATED  BELOW  3 

C 

PRINT  51 

51  F0RMAT[///.44X,4HlEFT,5X,5HRIGHT»5X.5HRIGHT,6X,4HLEFT,/,43X, 

15Ht-R0NT,5X.5HFR0NT,6X,4HREAR,6X»4HREAR,//»4X,5HTHETA,6X,3HPHl,6X. 

25HSL0PE.15X.4HF ( 1 ) ,6X . 4HF I ? I . 6X . 4HF [ 3  I . 6X. 4HF [4] , /.4X.5H [ DEG  3 ,5X. 
35HIDEG3 .5X.5HIPCT] ,14X,5HILBS3.5X,5H[LBS3 .5X.5HtLBS).5X.5HlL8Sl . 
4//  ] 

C  TH  =  DOWNSLOPE  DIRECTION  .  DEG  •  PI  /  180. 

C  PHI  =  SLOPE  ANGLE  .  DEL  «  INITIAL  INCREMENT  IN  PHI 

C 

DEL  =  5.  *  DEG 
073  =  .  75*DEL 

C  VAFY  TH  FROM  0  TO  360  DEGREES 

HO  650  J  =  l»  361 
AJ  =  J 

TH  =  [AJ  -  1.  1  *DEG 
C  ROTATE  GROUND  PLANE 

CAlL  ROTATE 

c  initialize  phi 

PHI  =  DEL 

OPH  =  DEL 

C  TILT  GROUND  PLANE 

2 Of)  CALL  TUT 


non  oo  n  on  norm 


CALCULATE  THE  4  WHEEt  -LOADS  TOR  THE  PRESENT  ORIENTATION 
AND  SLOPE. 

CALL  FORCE 

DETERMINE  WHETHER  INSTABILITY  EXISTS 
no  30n  I  =  1,4 
Ft  =  Fill 

IFIFF]  400*400,300 
300  CONTINUE 

IFIDPH  -  . 75*DEL )  310*310*320 
310  npH  -  .5  *  DPH 
320  ph!=  PHI*  DPH 
'50  TO  200 

<MiG  IFi  -FF  -  i.E-2  1  500.450,450 
<*50  DPH  s  .5  *  DPH 
PHI  =  PHI  "  OPH 
GO  TO  200 

500  S  =  1 0  0 . *S I NF [ PH  I  ]  /  COSFIPHI) 

REHFM8ER  THAT  DEG  =  PI/180. 

THU  =  TH  /  DEG 
PHI)  =  PHI/  DEG 

PRINT  THE  ORIENTATION  ANGLE.  SLOPE  ANGLE.  SLOPE.  AnD  THE 
FOUR  WHEEL-LOADS. 

PRINT  400 ,  THD.PHD.S,  lFtn,m*4} 

600  FOkPAT[3F10.2,10X*4Fi 0.2} 

65'.  CONTINUE 

RE  I  URN 
ENU 
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SUBROUTINE  ROTATE 

ROTATES  X » Y-PLANE  THROUGH  ANGLE.  AND  FINDS  NEW  COMPONENTS 
CF  ALL  VECTORS  NEEDED  FOR  WHEEL-LOAD  CALCULATION 

X » Y » Z  *  OLD  COORDS.  XP.YP.ZP  «  NEW  COORDS 

SUBROUTINE  ROTATE 

COMMON  STORAGE  FOR  STABILITY  STUOY#  ME-010 

Common  AJT131 ,  ALi ( 3 1 •  AL2I31.  CF131,  CGI3].  CLD13I.  CR l 3 ) . CRAX t 33 
1.  CXF(3],  CXRI3J.  F(<),  FAX133.  FLI31.  FLPI33,  TRtSl.  FRPI3) 

COMMON  OJP 1 33 »  0 JT (33.  R[9,3I,  RL13I.  RLPI33.  PRI33.  RRPI31.  Ut3l. 
lx[9},  XPC91.  XPPI93,  Y l 9 1 »  YP191,  YPP [93.  ZI9).  ZPI9I.  ZPPi93 

C 

COMMON  ALPHA,  DEG,  OA,  OSC,  PHI,  SHIFT,  T,  TH,  TPH,  W,  HI.  W2,  WF , 
lWLU,  WR,  WRAX,  WXF,  WXR 

C 

COMMON  IREAD,  NOSC.  NSTEER 
C 
C 

00  100  1=1.9 
C 

XPlI]  =  X [ I  1 *COSF [ TH 3  *  YtI]*SlNFtTHl 
YH  [  I  J  s  -X [ I ] *S I NF  t  TH ]  *  Y  1 1 1 *COSF i TH ] 

mo  zMin  =  zin 
c 

RE i URN 
ENU 
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SUbkOUTINE  7iLT 

TILTS  XP-AXIS  THROUGH  ANGLE  PH!  TOWARD  2P-AXIS*AND  FINDS  NEW 
COMPONENTS  OF  ALL  VECTORS  NEEDED  for  hheel-load  calculations 
XP, YP.ZP  S  OLD  COOROS*  XPP * YPP, ZPP  *  NEW  COORDS 
SUBROUTINE  TILT 

COMMON  STORAGE  FOR  STABILITY  STUDY,  ME-010 

COMMON  AJT13I.  ALU31.  AL2I31,  CFI31.  CGI3),  CLD13I#  CRl 3) ,CR4Xl3) 
i.  CXFI3I,  CXR (31*  F 1 4 J ,  FAXtSJ,  FL13].  FLPI3J,  FR[3J,  FRPI3) 

COMMON  0JP131.  OJTI31,  R l 9, 31 »  RLI3I,  RLPI31,  RRI3).  RRPI33,  UI31. 
1 X  I V  J ,  xP  1 9 ) »  XP?  19),  Y [ 9 1  *  YPI9I,  YPPI9 1 ,  ZI9),  JPI93,  ZPP19) 

COMMON  ALPHA,  DEG,  OA,  OSC,  PHI,  SHIFT.  T,  TH,  TPH,  W,  Hi,  W2,  WF. 
1HLU,  HR,  WRAX.  KXF .  WXR 

COMMON  IREAU,  NOSC,  NSTEER 


DO  100  1*1.9 
C 

XPKII]  *  XPH  1*C0SF  tPHl)  ♦  ZPin*SlNFlPHl) 
YPP  {  I  ]  *  YPUJ 

1  0  tt  7P?  [  I  1  *  -XPtD^SINF  (PHI)  ♦  ZPII  1*C0SFIPHI  } 
C 

RE  I  URN 
ENu 


>  2 
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SUHR0UTIN6  r0RCE 

CALCULATION  OF  THE  4  WHEEL  LOADS 
SUdROUTlNE  FORCE 

COMMON  STORAGE  FOR  STABILITY  STUDY,  HE-OlQ 

COMMON  AJTI3J.  AL1 ( 3 ] ,  AL2I3).  CFt3).  CG13),  CLD 1 3 5 •  CP [ 3 1 , CRAX ( 3 3 
1.  CxF  1 3  3  »  CXRI33,  FI4),  ^Axt3)»  FU31.  FLP(31,  FR{33,  FRPI31 
COMMON  OJP i 3 1 ,  0JTI31,  R{9,3),  RL 1 3  J  *  RLP  1 3 1 ,  RR;31»  RRP 1 3 1 »  U  t  3  J . 
IX { 9 ] ,  x?  f  9  3 ,  XPPI91,  Y  1 9  J «  YP 1 93 »  YPPI9),  ZI9),  ZPI93,  7PPt91 

COMMON  ALPHA,  DEG,  OA,  OSC.  PHI,  SHIFT,  T,  Th,  T*»H,  H,  HI,  W?,  hf, 
1HLO,  HR,  HrtAX,  WXF,  WXR 

COMMON  IREAD,  NOSC,  KSTEER 

DIMENSION  0141 

COMPUTE  COEFFICIENTS  OF  THE  FIRST  EQUATION 
01  =  XPPI9J*[YPPtlJ-YPPl8n  -  YPPI93*IxPPtl3-XPPt83 1 
02  =  XPP19J*IYPPI21-YPP[81 3  -  YPP [ 9 ] * [ XPP I  2 3 -XPP 1 8 3 3 
DC  1 =  XPP 1 9  J  * l YPP 1 5  3 -YpP  1  8  3  3  -  YPP i 9 3  *  I XPP 1 5 3  * XPP 1 8 3 3 

compute  determinant  of  coefficients 

OD  =  Dl*t lXPPl33-XPP[2Il*tYPPt43-YPPl23 1  -  I XPP t 4 3 -XPP  1 2 3)  * 

1  I YpP l 3 3 -YPP [  2 3  3  3  -  D2*l  t XPP 1 3 3 -XPP 1 1 31 • I YPP f 4 3 -YPP 1 1 ] ] 

2  -  tXPP(4l-XPPtl3 )*IYP°I3)-YPP’1)  3  3 

compute  denominator  determinants 

om  =  Wl*DCl*l  IXPPI3I-XPP121 3*tYPPt43-YPP[23  J  -  [  XPP  I  *  ) -XPP  f  ?  3  3  * 

1  l T  Pp I  3 ] -YPP [211  3  -  W*D2*l  lXPPI31-XPP[7j3*tYPPI43-YPPt7n  - 

2  UPPU1-XPP171  3*IYPPI31-YPPI713  3 

Of?]  =  W*D1*I  IXPP133-XPP  1  7  3  3  *[YPP[41-YPPI  7  3  3  -  I XPP l 4 ] -XPP l 7  ]  i  ♦ 

1  [TPPI3J-YPP171 33-  H1*DC1*[  IXPPI33-XPP11J 1  *  I YPP [ 4 1 -YP? [ 1 } 3  - 

2  [XPP[4)-XPPlll 1 •[ YPP I  3  3 -YPP { 1 3 3  3 

HU]  =  W*D1*  I  IXPP 1 7  3  -  XPP  12)1  *  t  YPPI 4  3  -  YPP  12  3 1  -  I  XPP  1 4  ] -XPP  I  2  3  3  * 

1  f YPP173 -YPPI2) J  3  -  W*D2*  l  iXPP[71-XPP(ll 3  * l YPP I  4 1 -YPP 1 1 3 :  - 

2  IX»P[4]-xPPm  J*[YPPf7)-YPPtl)  3  ]  ♦  Wl*DCl*  l  I  XPP  1 2  3 -XPP  1 1 3  3  * 

3  I YPP I  4  J - YPP 1 1 J  3  -  lXPPt4)-XPPtl3)*tYPPl21-YPPtl) 3  3 

0(4]  =  w*Dl*  I  IXPPI 33 -XPP123 )*[YPPI7J-YPP12) 1  -  I XPP l 7 3 -XPP r 2 1 1  * 

1  [YPPf 3 3 -YPP 1213  )  -  W*D2*  l  IXPP l 3 3 -XPP 111  1  * l YPP 1 7 ) -YPP { 1 33  - 

2  l  xpPi7)-xPPtl)  3*1  YPPI33 -YPPUll  3  *  Wl*DCl*  I  [  XPP  1 2  3 -XPP  1 1 3  3  «• 

3  f YpP!33-YPPf 11 3  -  [XPP133-XPPI1J 3*[YPP121-YPPI1)  )  1 

CALCULATE  THE  FOUR  WHEEL  LOADS 

no  100  I=l»4 
IOC  f'(  i  1  =  I)  I  I  J/DD 
C 


RE l URN 
END 
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C  HDOTVR 

SUHROUTINE  MDOTVR 

THIS  SU3R0UT f  ‘>iE  COMPUTES  CRITICAL  SLOPES  AND  ORIENTATIONS  FOR 
TURNOVER  FOR  A  VI;  ICLE  WITm  MIORANGE  OSCILLATION.  IT  IS  ASSUMED 
THAT  THE  VEHICLE  IS  INITIALLY  SITTING  ON  A  (HORIZONTAL!  PLANE  IN 
A  COMPLETELY  OSCILLATED  CONDITION  WITH  ONE  OF  THE  A  TIRES  OFF  THE 
PLANE  AND  THE  OTHER  3  ON  THE  PLANE.  TURNOVER  CALCULATIONS  ARE 
HAUF  FOR  each  of  the  4  TIRES  OFF  THE  PLANE. 


COMMON  STORAGE  FOR  STABILITY  STUDY,  ME-Q10 

COMMON  A JT 1 3 ] ,  ALl [ 3 ] ,  Al2(3!,  CFI3I,  CG(3l.  CLDI3I,  CR 1 3 ) , CRAX 1 3 J 
1.  CXFI3I,  CXR { 3 ! *  F 1 4 1 ,  F  AX  I  3 ! »  FLI31,  FlP[35»  FRI3J ,  FRPI31 
COMMON  0JPI3J,  0 JT  13),  R I  9, 3 1 ,  RL131.  RLP(3J,  RRI3I,  RPPI3),  U tc ! , 
IX l v ) ,  XP [ 9 i ,  XPP (91,  Y  t  9 ) ,  YP I  9 ) ,  YPPt9l,  Z(9),  ZP(9),  ZPPI9) 

COMMON  ALPHA,  deg,  GA,  OSC,  PHI,  SHIFT,  T,  TH,  TPH,  H,  HI*  W2,  WF , 
1HLU,  HR,  HR AX ,  WXF,  WXR 

COMMON  iread,  NOSC,  nsteer 

"I  MENS  I  ON  Cl(3],C2(3],CCl(3) , CC2I3J ,FC(3! 

PRINT  1 

1  format i iHi, 2ox , t>iH turnover  assuming  complete  oscillation  has  occur 

1REU  ] 

TAU  =  -OSC  *  DEG 

00  mo  i si , 3 

Cl l I  1  a  P(5,I! 

100  C2UJ  S  R  [  6,  I  ] 

DO  7000  I  =  1,4 

DEL  =  !-l.J**(I*l) 

30  TO  (1000, ?000, 3000, «000l  ,  I 

TURNOVER  TO  LEFT,  KITH  RIGHT  FRONT  TIRE  OFF  GROUND 

100<>  CALL  SROTIOJP.U,!  AU,FLP*FC) 

CALL  SROT (0JP.U,TAU,C1 .CC1) 

C 

CALL  OSNlRHP,RLP.*FC,DEL,RLP,  ALl! 

CALL  OSN(RRP,RLP.FC,DEL»FC 
CALL  OSN(RRP,RLP,FC.DEL,CC? ,CC1) 

CALI  OSN(RRP,RLP»FC,DEL,C2  ,CC2! 

C 

PRINT  ? 

2  FORMAT r///,10X,54H(A!  TURNOVER  TO  LEFT,  WITH  RIGHT  FRONT  TIRE  OFF 
1GROUND  ,///] 
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GO  TO  SHOO  __ rtl ...n 

;  TURNOVER  TO  LEFT ,  WITH  RIGHT  REAR  TIRE  OFF  GROUND 

2000  CALL  SROT!OJP,U.TAU#RLP»FC1 
call  SR0TI0JP.tI.TAU.C2.CC2] 

CALL  OSMFRP.TLP.FC.nEL.FC  »AL1) 

CALt  OSMFRP.FLP.FC.P6L.FLP.AL2) 

CALL  OSMFRP.FLP.FC.'JfcL.Cl.CCl  J 
.j  Ai_L  OSUFRP.FLP.FC.PEL.CC2.CC2] 

C 

i  roRMAT^///.10X.5SHtb]  TURNOVER  TO  LEFT.  WITH  RIGHT  REAR  TIRE  OFF  G 
1M0UN0  .///I 

r* 

O 

c  ';0  TO  RIGHT, UITH  LEFT  FRONT  TIRE  OFF  GROUND 

c 

TA)  =  -TAU 

s- alL  SROTIOJP.tJ.TAU.FRP.FCl 
CAU  SROTIOJP.U.TAU.Ci.CCll 

£ 

CALL  OSNlRRP.RLP.FC. DEL. RRP.ALU 
CAU  OSNlRHP.RLP.FC.PEL.FC  »AL21 
■'all  OSN  l  RRP » RLP » FC .  *)EL » CCl »  CC1 1 
CALL  OSN l RHP » RLP » FC . PEL  * C2  *CC2) 

C 

-  rOKhIn///.10X.5«H[C)  TURNOVER  TO  RIGHT.  WITH  LEFT  FRONT  TIRE  OFF 
1CR0UNP  ,///) 

C 

'50  Turnover  to  right,  -ith  left  re.r  tire  off  ground 

Atiui.  CALI  SROTIOJP.U.TAU.PRP.FC  j 
CALL  SHOT ( OJP »U » T AU* C2  »CC2l 

v:ai.L  OSH  l  fkp # F  LP » F  C ♦  PEL » r  C,  AL1  I 
■"'ALi.  0SNlFHP.FLP.FC.DEL.FRP.AL2) 

CALL  OSNiFRP.FLP.FC»PEL»Cl.CCl  ] 

'.AlL  0SNtFKP.FLP.FC.PEL.CC2.CC2) 

C 

rSRNATU//.10X.S3H!DJ  TURNOVER  TO  RIGHT.  WJTK  LEFT  REAR  TIRE  OFF  G 

1R0UM)  ,///) 


CO-PUTfc  COMBINED  CG  FOR  OSCILLATED  VEHICLE 


Si; Li1''  00  SlOH  jsl.3 

*J  1 » [•  vGij)  =  (  CC1I  JI*W1  *  CC2 1 J )  *W2  J  /  W 
CAlI  SLOPE 

iFiPHl  -  1.6-4]  7000.5200*5200 
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c 

D r'. 0 i*  GO  TO  I5500#5600. 5500, 5600) # 1 
C 

t>5l)l.  IF  U-S7EER-2)  6  000,5550,6000 
C 

•5iJbO  OA  =  OA  -  ALPHA 
GO  TO  6000 
C 

56G0  IFiNSTEER-U  6000,5650,6000 
C 

G6t>n  OA  =  OA  ♦  ALPHA 
C 

6 U 0 1)  PRINT  6,  PH»,TPH,OA,T  _  „  ^ 

a  rOKMATI/-F20.?»6H  DEG  =»F9.2»4H  PCT»F21-2»4H  DEG# F51 . 5 . 2H=T » //// 1 

C 

7001*  CONTINUE 
C 

re  turn 
fn  J 
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APPENDIX  V 


•*••»«<«%*  ««  /vwmtirrtct  t»m  f* ACP  T  AMI; 
rj\nmrLE.  uuiruio  av«i  ■*  *“■“ 


ii  yranrn  'is 


Following  is  the  output  from  three  test  runs  of  the  computer  program. 
The  first  test  run  is  a  stability  study  •  '  a  two-inch  cube  weighing  ten 
pounds.  (The  computer  assumes  the  cube  a  vehicle  with  the  wheel  loads 
located  at  the  bottom  four  comers  of  the  cube. ) 


The  second  and  third  test  runs  are  stability  studies  of  a  fictitious 
vehicle  with  rear  axle  oscillation.  The  second  run  uses  the  action  line 
method  and  the  third  run  uses  the  wheel  load  method.  In  the  third  run,  the 
computer  assumes  that  the  vehicle  has  midrange  oscillation  with  the 
oscillation  joint  located  at  the  center  of  the  rear  axle. 


TEST  RUN  1  •  static  STABILITY  STUDY  OF  A  CUBE 


vehicle  data 


l  2.0,  1.0*  0.0]  «rL 

t  2.0,  -1.0#  0.0)  «FR 

t  0.0,  -1.0#  0.0)  «RR 

I  0.0,  1.0#  0.01  »RL 

[  0.0,  0.0#  0*01  »AJT 

t  0.0,  0.0#  0.0)  »DJT 

(  0.0,  0.0#  0.0)  »FAX 

I  r.0,  0.0.  1'0)  «CF 

(  0.0,  0.0#  0.0)  aCLD 

l  0.0,  0.0#  0.0)  »CXF 

j  0.0,  0.0#  0*0)  *CR 

{  0,0,  0.0#  0.0)  *CRAX 

t  0.0#  0.0#  0*0)  •  »CXR 

0.00  sALHHA  (DEG)  ■  ANGLE  OF  LEFT  STEER 
0.00  bOSCIDEG) >MAX  OSCILLATION  ANGLE 
0.00  *  SHIFT#  MEASURED  FROM  CENTER  OF  TREAD 
3  aSTEER  TYPE11*ARTICULATED»2«REAR  WAGON# 3*ACKERNANN) 

1  »0SCILLATI0N  TYREU«N0  0SC»2«REAR  AXLE  0SC#3*MlDRAN0E  080 


10.0  swF 
0.0  =WLC 

o.o  *wxr 

0.0  bWR 
0.0  sWRAX 
0.0  sWXR 
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OOWHSLOPE  tS  TO  THE  LEFT 


TEST  RUN  24  -  REAR  AXLE  OSCILLATION  -  USES  ACTION  LINE  METHOD 


VEHICLE 

DATA 

100.0, 

30.0, 

0.0) 

aFL 

100.0, 

-30.0, 

0.0) 

■  PR 

0.0, 

—36.0, 

0.0) 

■  RR 

0.0, 

30.0, 

0.0) 

■RL 

50.0, 

0.0, 

30.0) 

■  A  JT 

l  O-.C,  0,0.  30.8)  'OJT 

t  100.0.  0.0.  30.0)  aPAX 

(  6o.o,  o.o,  40.0)  ter  sooo.o  «wr 

(  150.0,  0.0,  50.0)  aeLD  2000.0  *WLD 

t  0.0,  0.0-  0.0)  afixr  0.0  *MXr 

(  0.0,  0.0,  0.0)  BCR  0.0  aWR 

(  0.0,  0.0,  30.0)  aCRAX  1000.0  «WRAX 

t  O.o,  0.0,  0.0)  aCXR  0.0  aUXR 

0.00  'ALPHA  IDEQ1  ■  ANGLE '  OP  v.-fFT  STEER 
15.00  aOSC(DEG) «MAX  OSCILLATION  ANGLE 
0.00  •  SHIFT,  MEASURED  FROM  CENTER  OF  TREAD 

1  ‘STEER  TYPE t la ART  I CUL ATED» 2aREAR  WAGON, S'ACKERMANN I 

2  'OSCILLATION  TYPEll'NO  0SC,2aREAR  AXLE  0SC,3*MIDRANQE  080 


SIDFSLO’E  TURNOVER  ASSUMING  NO  OSC ! LI  AT ! ON 


S1DESLOPF  TURNOVER  ASSUMING  NO  OSCILLATION 


o 


TEST  RUN  20  *  H10RANGE  OSC  AT  REAR  AXLE  *  USES  WMfctL  lC»u  wETriGu 


vehicle  Data 


3  «OSCILLATION  TYR£(l«NO  OSC,?»REiR  AXLE  OSC,S*HlDRANBE  OSC) 


( 

85.7.  0 

.0.  42.9] 

«  FRONT  CG 

7000 

.1)  >  FRONT  WEIGHT 

t 

8.0.  G 

.0.  30.0) 

•  REAR  ce 

1000 

.0  •  REAR  WEIGHT 

t 

75.0.  0 

.0.  41.2) 

«  COMBINED 

CO  8000 

.0  «  TOT/.L  weight 

FOLLOWING  IS  A  TABLE  GIVING  SlORES  AT  WHICH  INSTABILITY  OCCURS  FOR  A  VEHICLE  1 

WITH  midrange  oscillation 

THETA 

*  DOWMSLORE 

DIRECTION  (MEASURED  COUNTERCLOCKWISE  FROM  I 

ROS  X-AXtS) 

RHI 

*  ANGLE  THNU 

WHICH  GROUND 

RLANE  must 

BE  TILTED  FOR  INSTABILITY  TO  OCCU* 

SLOPE 

*  T  AN ! RH I )  I 

>:  rercent 

Fill 

*  WHEEL  LOAD  AT  TILT  ANGLE  RHI,  AS 

indicated 

below 

left 

RIGHT  RIGHT 

Left 

front 

FRONT 

REAR 

REAR 

THETA 

PHI 

SLOPE 

F 1 1 1 

r  1 2 1 

F 1 3 ) 

r  1 4  j 

tDEG] 

[DEG] 

[PCT1 

(LBS) 

(LBS)  (LBS) 

(LBS) 

0.00 

31.22 

60.61 

4000.01 

4000.01 

-0.01 

-0.01 

1.00 

31.22 

60.62 

4058.19 

3941.82 

-0.00 

-0.00 

2.00 

31.23 

60.64 

4116.41 

3883.61 

-0.01 

-O.oi 

3.00 

31.25 

60.69 

4174.70 

3825.31 

-0.00 

-0.00 

4.00 

31.28 

60.75 

4233.10 

3766.91 

-0.01 

-0.01 

5.00 

31.32 

60.84 

4291.64 

3708.37 

-o.oi 

-0.01 

6.00 

31.36 

60.94 

4350.36 

3849.66 

-0.01 

-o.oi 

7.00 

31.41 

61.06 

4409.30 

3590.72 

-0.01 

-0.01 

8.00 

31.47 

61.20 

4468.49 

3531.54 

•0 . 01 

-0.01 

9.00 

31.53 

61.36 

4527.96 

3472-05 

•0.01 

-o.oi 

10.00 

31.61 

61.54 

4587.76 

3412.24 

-0.00 

-0.00 

11.00 

31.69 

61.74 

4647.95 

3352.07 

-0 . 01 

-0.01 

12.00 

31.78 

61.96 

4708.53 

3291.48 

-8.01 

-0.01 

13.00 

31.88 

62.20 

4769.57 

3230.44 

-0.01 

-0.01 

14.00 

31.99 

62.46 

4831.11 

3168.91 

-6.01 

-o.oi 

13.00 

32.11 

62.74 

4893.18 

3106-84 

-o.oi 

-0.01 

16.00 

32.23 

63.05 

4955.82 

3044.18 

-b.oo 

-0.00 

17.00 

32.36 

63.38 

5019.12 

2980.90 

-8.01 

“0.01 

18.00 

32.51 

63.72 

5083.07 

2916.93 

-0.00 

-0.00 

19.00 

32.66 

64.10 

5147.77 

2852.24 

-0.00 

-0.01 

20.00 

32.82 

64.50 

5213.25 

2786.76 

-0.01 

-o.oi 

21.00 

32.99 

64.92 

5279.57 

2720-45 

-0.01 

-0.01 

22.00 

33.17 

65.37 

5346.77 

2653.24 

-0.00 

-0.01 

23.00 

33.36 

65.84 

5414.95 

2585.00 

-0.01 

-0.01 

24.00 

33.56 

66.34 

5404.11 

2515.90 

-0.00 

-0.01 

29.00 

33.77 

66.87 

5554.30 

2445.64 

-0.00 

-0.01 

26.00 

33.99 

67.43 

5625.01 

2374.21 

-e.oi 

-0.02 

27.00 

34.22 

68.02 

5690.44 

2301*57 

-0.00 

-0.01 

28.00 

34.47 

68.64 

5772.36 

2227.64 

-0.00 

-0.00 

29.00 

34.72 

69.29 

5047.71 

2152.30 

-0.00 

-0.01 

30.00 

34.99 

69.98 

5924.54 

2(575.40 

-0.01 

-0.02 

31.00 

35.26 

70.71 

6002.09 

1997.12 

-0.00 

-0.01 

32.00 

35.55 

71.47 

6002.91 

1917.10 

-*.00 

-0.01 

33.00 

35.85 

72.27 

6164.74 

1835.29 

-0.01 

-0.02 

34.00 

36.17 

73.10 

6240.39 

1751-03 

-0.00 

-0.01 

39.00 

36.50 

73.99 

6334.05 

1665.96 

-0.00 

-0.01 

36.00 

36.84 

74.91 

6421.05 

1570.17 

-0.01 

-0.02 

37.00 

37.19 

75.89 

6511.91 

1408.12 

-0.01 

-0.02 

38 .00 

37.56 

76.91 

6604.30 

1395.70 

-0.00 

-0.01 

39.00 

37.95 

77.99 

6699.31 

13a0*71  _ 

-o.oo 

-0.01 

40.00 

38.35 

79.12 

6797.02 

1202.99 

-0.00 

-0.01 

41  .00 

38.77 

80.30 

6097.63 

1102.37 

-0.00 

-0.00 

42.00 

39.20 

81.55 

7001.37 

990.64 

-0.00 

-0.01 

43.00 

39.65 

82.87 

7100.46 

091.57 

-0 » 01 

-0.03 

44.00 

40.12 

84.25 

7219.04 

780.99 

-0.01 

-0.02 

49.00 

40.60 

85.71 

7333.36 

666.65 

-0.00 

-0.01 

46.00 

41.10 

87.25 

7451.00 

540.21 

-0.00 

-0.01 

47.00 

41.63 

88.87 

7574.67 

425.37 

-0.01 

-0.03 

48.00 

42.17 

90.58 

7702.17 

297.08 

-0.01 

-0.04 

49.00 

42.73 

92.38 

7034.66 

165.30 

-0.01 

-0.03 

50.00 

43.32 

94.29 

7972.62 

27.42 

-0.01 

-0.03 

51.00 

43.04 

93.39 

7939.45 

-0.00 

8.30 

52.25 

52.00 

42.57 

91.87 

7066.46 

-0.01 

10.44 

115.12 

53.00 

42.12 

90.42 

7795.71 

-0-01 

28.39 

175.91 

54 .00 

41.68 

89.04 

7727.07 

-0-00 

30.17 

234.77 

55.00 

41.26 

07.72 

7660.45 

-o.oo 

47.70 

291.79 

56 .00 

40.85 

86.47 

7595.66 

-0.00 

57.24 

347.10 

57.00 

40.46 

85.27 

7532.65 

-0.00 

66.56 

400.00 

58.00 

40,08 

84.13 

7471.30 

-0.01 

75.74 

452.97 

59.00 

39.71 

83.05 

7411.51 

-0.00 

04.79 

503.71 

60 .00 

39.36 

02.01 

7353.10 

-0.00 

93.71 

553.11 

61  .00 

39.02 

81.02 

7296.26 

-0.01 

102.53 

601.23 

62.00 

38. ^ 

80.08 

7240.62 

-0.00 

111.23 

640.15 

63.00 

38.37 

79.18 

7106.23 

-0.01 

119.83 

693.94 

64.00 

38.07 

78.32 

7133.00 

-0.00 

120.34 

738.66 

69.00 

37.78 

77.50 

7000 . 80 

-0.00 

136.76 

782.36 

66.00 

37.50 

76.72 

7029.79 

-0.00 

145.10 

025.11 

67.00 

37.23 

75.98 

6979.60 

-0.00 

153.36 

066.96 

68.00 

36.97 

75.27 

6930.50 

-0.00 

161.55 

907.95 

69.00 

36.72 

74.60 

6002.20 

•0.00 

169.67 

948.13 

70.00 

36.49 

73.96 

6834.72 

-0.00 

177.73 

987.55 

71.00 

36.26 

73.35 

6788.02 

-0.00 

105.73 

1026.25 

72.00 

36.04 

72.77 

6742.05 

-0.00 

193.69 

1064.26 

73.00 

35.84 

72.22 

6696.78 

-0.00 

201.60 

1101.63 

74.00 

35.64 

71.70 

6652.15 

-0.00 

209.46 

1130.39 

79.00 

35.45 

71.20 

6608.15 

-0.00 

217.29 

1174.56 

76.00 

35.27 

70.74 

6564.72 

-0.00 

225.09 

1210.19 

77.00 

35,11 

70.30 

6521.83 

-0.00 

232.06 

1245.31 

7S  =  00 

79.00 


Si.  0,0 
82.0,0 
83 . 0,0 
84.00 
85*00 
86.00 
87.00 

68.58 
89.00 
90.00 
91.00 
92.00 
93.00 
94.00 
95.00 
96.00 
97.00 
98. 00 
99.00 
100.00 
101.00 
102.00 
103.00 
104.00 
105.00 
106.00 
107.30 
108.00 
109.00 
110.00 
111.00 
112.00 
113.00 
114.00 
115.00 
116.00 
117.00 
118.00 
119.00 
120.00 
121.00 
122.00 
123.00 
124.00 
125.00 
126.00 
127.00 
128.00 
129.00 
130.00 
131.00 
132.00 
133.00 
134.00 


34.9? 

34.79 

34.6? 

34.52 
34.39 

34.28 

34.17 
34.07 

33.98 

33.89 

33.82 

33.75 

33.69 

33.64 

33.59 
33.56 

33.53 

33.50 

33.49 

33.48 

33.48 

33.49 

33.51 

33.53 
33.56 

33.60 

33.64 

33.70 

33.76 

33.83 

33.90 

33.99 
34.08 

34.18 

34.29 

34.41 

34.53 
34.66 

34.80 
34.95 
35.11 

35.28 

35.45 

35.64 

35.83 
36,03 

36,24 

36.46 

36.69 
36.93 

37.18 
37.44 

37.71 

37.99 

38.29 
38.59 

38.90 


69.88 

69.49 
69.12 

68.76 

68.46 

66.16 

67.88 
67.63 
67.40 

67.18 

66.99 
66.82 
66.67 

66.53 

66.42 

66.33 

66.25 

66.20 

66.16 

66.14 

66.15 

66,17 

66.21 

66.26 

66.34 
66.44 
66.55 

66.69 

66.84 
67.01 
67.21 

67.42 
67.66 

67.91 

68.19 

66.48 

68.80 
69.14 

69.51 

69.90 
70.31 

70.74 

71.20 

71.69 

72.20 

72.74 
73.30 

73.90 

74.52 

75.18 
75.86 
76.58 
77,33 
7f  .11 
78.94 
79.79 

80.69 


6479.46 

-n .  01 

6437.55 

•  0  •  00 

6396*10 

-0-00 

6355.06 

-0.03 

6314.41 

-0.01 

6274.13 

-0.01 

6234.17 

-0-01 

6194.52 

-0.00 

6155.15 

-0.01 

6116.03 

-o.oo 

6077.16 

-0.01 

6038.49 

•0.00 

6000.00 

-0.00 

5961.68 

•  o.oo 

5923.51 

•0.00 

5885.45 

-0.01 

5847.49 

-0.31 

5809.61 

-0.01 

5771.78 

-0*00 

5734.00 

-0.01 

5696.21 

-0-00 

5656.44 

-0-01 

5620.62 

-o.oo 

5582.76 

-3.01 

5544.83 

-0.03 

5506.81 

-2.00 

5468.69 

-0.01 

5436 . 43 

-0^00 

5392.02 

-0.01 

5353.43 

-o.oo 

5314.64 

-  0 . 1*  1 

5275.64 

-0.01 

5236.39 

-0.01 

5196.87 

-o.oo 

5157.07 

•  o.oo 

5116.95 

-o.oo 

5076.50 

-0.01 

5035.68 

-0.00 

4994.47 

-o.oo 

<952.85 

-0.01 

4910.77 

-3.00 

4868.23 

-0.01 

4825.18 

-0.00 

4781.60 

-0-00 

4737,45 

-0-00 

4692.70 

-o.oo 

4647.32 

-o.oo 

4601.27 

-o.oo 

4594.53 

-0.01 

4507.03 

-o.oo 

4458.76 

-o.oo 

4409.66 

-0*00 

4359.70 

-0.00 

4308.83 

-o.oo 

4257.00 

-o.oo 

4204.17 

-0.00 

4150.28 

-0.30 

240.61 

1279.94 

24e . 33 

1314.12 

256.05 

1347.86 

263 . 75 

1551 . 19 

271.45 

1414.14 

279.14 

1446.74 

286.84 

1479.00 

294.55 

1510.94 

302.26 

1542.59 

310.00 

1573.97 

317.75 

1605.10 

325.53 

1635.99 

333.33 

1666.67 

341.17 

1697.15 

349.05 

1727.45 

354.97 

1757.59 

364.93 

1787.56 

372.95 

1817.44 

381.03 

1847.19 

389.17 

1876.84 

397.38 

1906.41 

405.66 

1935.92 

414.02 

1965.37 

422.46 

1994.79 

430.99 

2024.16 

439.62 

2053.57 

448.35 

2082.96 

457.20 

2112.38 

466.15 

2141.84 

475.24 

2171.34 

484.45 

2200.91 

495.80 

2230.57 

503.30 

2260.31 

512.96 

2290.17 

522.78 

2320.15 

532.78 

2350.27 

542.96 

2380.54 

553.34 

2410.98 

56"  93 

2441.61 

574 ,  /  o 

2472.43 

565.76 

2503.47 

597.04 

2534.74 

608.58 

2566. 23 

620.39 

2598.02 

632.49 

2630.07 

644.89 

2662.41 

657.62 

2695.06 

670.69 

2728.04 

684,12 

2761.36 

697.94 

2795.03 

712.16 

2829.08 

726.82 

2863.52 

741.94 

2898,37 

757.54 

2933.63 

773.67 

2969.34 

790.34 

3005.50 

807.60 

3042,12 

6  h 


t 


135.0  0 
134.00 
137.00 
138.00 
139.00 
1*0.00 
141.00 
142.00 
143.00 
144.00 
148.00 
1*6. 0*’- 
147.00 
148.00 
149.00 
150.00 
151.00 
152-00 
153.00 
154.00 
155. 00 
158.00 
157.00 
158.00 
159.00 
160.00 
161.00 
162.00 
163.00 
164.00 
165.00 
166.00 
167,00 
168.00 
169.00 
170.00 
171,00 
172.00 
173.00 
174.00 
175.00 
176.00 
177.00 
178.00 
179.00 
180.00 
181.00 
182.00 
183.00 
184.00 
185.00 
186.00 
187.00 
188.00 
189.00 
190.00 
191.00 


39.22 

8i  .63 

39.56 

62.60 

39.90 

83.62 

40.26 

84.69 

40.63 

85.80 

41.01 

86.96 

41.40 

88.16 

4i.8G 

69.42 

42.22 

90.73 

42.65 

92.10 

43.08 

93.53 

43.53 

95.01 

44.00 

96.56 

44.47 

98.17 

44.96 

99.84 

45.45 

101.59 

45.96 

103.40 

46.48 

105.29 

47.00 

16“  .25 

47.54 

109.28 

48.08 

111.39 

48.64 

113.58 

49.20 

115.85 

49.77 

118.19 

50.34 

120.62 

50.92 

? 23.12 

51.50 

125.69 

52.08 

128.35 

bi .  66 

lui .07 

53.24 

133.85 

53.81 

136.70 

54.39 
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